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1.  IN  TRODUCTICN 


This  report  describes  the  development  and  use  of 
certain  computer  programs  which  allow  the  determination  of  rele¬ 
vant  information  with  respect  to  strength  calculations  for  all 
types  of  involute  gears  including  fuze  gears. 

Both  programs,  LKWIS  CIRCLES  and  LEWIS  ENVELOPE, furni sh 
the  coordinates  of  the  left  side  of  any  involute  tooth  (See  Figure  1) , 
together  with  the  associated  Lewis  and  AGMA  geometry  factors.  The 
determination  of  these  factors  is  based  on  the  assumotion  that  the 
maximum  contact  force  F  is  applied  at  the  tip  of  the  tooth  and  that 
therefore  the  origin  V  of  the  inscribed  Lewis  parabola  is  located  at 
the  point  of  intersection  of  the  line  of  action  of  this  force  and  the 
tooth  centerline.  The  point  of  tangency  between  this  parabola  and  the 
tooth  profile,  together  with  the  x-coordinate  XTAN  required  for  the 
determination  of  the  above  factors,  is  found  by  an  iterative  process. 

The  involute  portion  of  the  tooth  profile  is  obtained  in 
both  programs  by  a  mathematical  description  of  the  unwinding  of  the 
base  circle  tangent.  The  trochoid  or  root  portion  of  the  tooth 
profile  is  obtained  by  two  different  methods. 

Program  LEWIS  CIRCLES  finds  the  associated  coordinates 
by  a  simulation  of  the  hob  cutting  or  drafting  process,  while  program 
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LEWIS  ENVELOPE  makes  use  of  an  analytical  method,  based  on 
the  relationship  between  generating  curves  and  envelopes. 

While  the  results  of  both  programs  are  identical,  program 
LEWIS  CIRCLES  cannot  be  used  when  the  hob  fillet  radius  is 
zero . 

To  assist  the  user  as  much  as  possible  with  the  details 
of  operation  of  these  programs,  the  following  outline  of  input 
and  output  parameters,  is  given:  (The  input  and  output  parameters 
of  both  programs  are  identical) . 

a.  Input  Data  (All  given  for  unity  diametral  pitch) 

1 .  Number  of  Teeth 

2.  Effective  distance  (see  Appendix  B) 

3.  Unity  diametrical  pitch 

4.  Hob  pressure  angle 

5.  Circular  tooth  thickness  at  standard 

pitch  circle 

6.  Hob  fillet  radius 

7.  Gear  tooth  addendum 

8.  Angular  increments  for  involute  and 

trochoid  generation 


b.  Output  Data 


1.  Results  of  various  intermediate  computations 
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2. 


Lewis  factor 


3.  AGMA  geometry  factor 

4.  Various  parameters  associated  with 

the  point  of  tangency  between  Lewis 
parabola  and  tooth  profile 

5.  Coordinates  of  Tooth  profile 

As  a  further  aid  to  the  user,  material  on  which  this 
report  is  based  is  compiled  in  the  appendixes. 

Appendix  A:  Determination  of  Involute  Coordinates 

This  appendix  gives  the  derivation  of  the  x  and  y-coor- 
dinates  of  the  involute  portion  of  the  tooth  profile  of  an  arbi¬ 
trary  tooth  in  terms  of  the  coordinate  system  shown  in  Figure  1. 
Use  is  made  of  the  angle  a  which  describes  the  unwinding  of  the 
base  circle  tangent. 
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Appendix  B:  Determination  of  Path  Coordinates  of  Effective 

Rack  Point 


After  defining  the  concepts  of  effective  rack  point 
and  effective  distance  for  hobs  with  and  without  fillet  radii,  this 
appendix  gives  derivations  for  the  path  coordinates  of  the  effective 
rack  point  for  undercut  and  non-undercut  gear  teeth.  In  order  to 
accomplish  this,  the  generating  angle  i|j,  which  describes  the  rolling 
of  the  rack  pitch  line  on  the  gear  pitch  circle,  is  introduced.  Again, 
the  results  are  given  in  terms  of  the  tooth  centerline  coordinates 
shown  in  Figure  1. 

Appendix  C:  Generation  of  Tooth  Root  in  Presence  of  Fillet 

Ra dius  on  Rack  Cutter 

Two  methods  of  obtaining  the  coordinates  of  the  tooth 
root  are  shown.  The  first  method  simulates  the  cutting  or  drafting 
process  and  is  used  in  program  LEWIS  CIRCLES.  The  second  method,  which 
is  based  on  certain  relationships  between  generating  curves  and 
envelopes,  forms  the  basis  for  the  root  geneiation  in  program  LEWIS 
ENVELOPE. 
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Appendix  P:  Initial  and  Final  Angles  for  Trochoid 

and  Involute  Generation 


This  appendix  first  gives  a  review  of  existing  procedures 
for  the  determination  of  the  inner  form  radius  for  undercut  and  non¬ 
undercut  teeth.  Subsequently,  expressions  for  the  limits  of  the 
generating  angles  a  and  are  derived  for  both  types  of  teeth. 

Appendix  E:  Gear  Strength  Calculations.  Determination 

and  AGMA  Geometry  Factors 

Initially,  the  origin  of  the  Lewis  gear  strength 
formula  is  reviewed  and  expressions  for  the  determination  of  various 
associated  parameters  are  derived.  Subsequently,  the  iterative 
procedure  for  finding  the  point  of  tangency  between  the  Lewis  parabola 
and  the  tooth  profile,  which  is  used  both  in  program  LEWIS  CIRCLES 
and  in  program  LEWIS  ENVELOPE,  is  explained. 

Pl9RIam  MATRIX 

This  appendix  describes  and  lists  program  MATRIX,  which 
is  used  to  determine  the  size  of  the  arrays  in  the  two  computer 
programs.  it  also  contains  sample  outputs  of  program  MATRIX 
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which  pertain  to  the  sample  runs  of  both  LEWIS  programs. 

(See  Section  3.) 

Appendix  G:  Program  LEWIS  CIRCLES 

A  complete  listing  of  the  program  LEWIS  CIRCLES  is 
given,  together  with  outputs  of  the  two  sample  runs. 

Appendix  H:  Program  LEWIS  ENVELOPE 

This  appendix  contains  a  complete  listing  of  program 
LEWIS  ENVELOPE.  In  addition, the  results  of  the  two  sample  runs 
are  shown. 
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DESCRIPTION  OF  COMPUTER  PROGRAMS 


The  following  section  gives  descriptions  of  the 
essential  steps  of  the  programs  LEWIS  CIRCLES  and  LEWIS  ENVELOPE. 

Program  LEWIS  CIRCLES  is  listed,  together  with  sample  runs, 
in  Appendix  G.  Appendix  H  furnishes  the  same  for  program  LEWIS 
ENVELOPE.  The  specific  sample  runs  of  both  programs  are  discussed 
in  Section  3. 

The  input  and  output  parameters  of  both  programs  are 
identical  and  will  therefore  only  be  discussed  in  detail  in 
conjunction  with  program  LEWIS  CIRCLES. 

Since  the  vector  Y  and  the  matrix  X,  as  used  for  the 
determination  of  the  root  profile  in  LEWIS  CIRCLES  (see  also  Appen¬ 
dix  C-2 :  Iterative  Method  of  Obtaining  Root  Profile),  may  be 
very  large,  it  is  necessary  to  use  care  in  the  dimensioning  of 
these  arrays  as  well  as  in  the  request  for  computer  storage  space. 
Appendix  F  gives  program  MATRIX  which  allows  the  determination  of 
the  sizes  of  Y  and  X,  as  well  as  of  the  vectors  XI  and  Y1  which  are  used 
for  the  initial  storage  of  the  involute  coordinates. 

I.  Program  LEWIS  CIRCLES  (See  Appendix  G) 

a^ _ Dimensioning  of  Arrays 

(For  more  details  on  the  use  of  the  various  arrays 
see  description  of  computations  below) . 
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The  vectors  XI  and  Y1  are  used  for  the  initial  storage  of 
the  involute  coordinates  of  the  gear  tooth  (see  Appendix  A) . 

The  minimum  number  of  locations  for  either  of  these  arrays 

is  given  by:  (See  also  Appendix  F) 


1 1 


fin 


a . 
in 


Aa 


tt) 


where 

ou n  =  initial  involute  generating  angle  [See  equ.  (D23) 
of  Appendix  D] . 

=  final  involute  generating  angle  [See  equ.  (D2  4)]. 

\a  =  increment  of  involute  generating  angle. 

The  arrays  X2  and  Y2  are  used  for  the  final  storage  of  all 
values  cf  the  tooth  profile  coordinates.  The  required  number  of 
locations  is  represented  by  the  sum  of  those  needed  for  the 
vectors  Y1  (or  XI)  and  Y.  It  can  only  be  determined  once  the 
length  of  the  vector  Y  is  known.  The  later  is  obtained  with  the 
help  of  the  aforementioned  program  MATRIX,  which  also  finds  the 
size  of  matrix  X. 
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_ Input  Parameters 


The  following  parameters  represent  the  required 
input  data  for  the  program.  (They  are  also  printed  out  on 
top  of  the  first  output  page). 


N 

BEFF  = 

PD 

THETAD 


number  of  teeth  in  gear  under  consideration 
b^^  for  unity  diametral  pitch  (See  Appendix  B) 

=  1.000,  unity  diametral  pitch  of  hob 

0,  the  pressure  angle  of  the  hob  (the  D  at  the  end 
expresses  that  an  angle  is  given  in  degrees) 


CAPTCS  =  Tcs ,  the  circular  tooth  thickness  of  the  gear  at  the 

standard  pitch  radius  R  =  yp— ,  in  terms  of  Pd  =  1.000 

P  d 

RC  =  r  ,  the  rack  fillet  radius  for  unity  diametral  pitch 

KADD  =  Addendum  constant  of  gear,  with  the  magnitude  of  the 

addendum  defined  by  KADD/PD  ,  in  terms  of  Pd  =  1.000 


DELPSI  =  Aip,  the  increment  of  the  roll  angle  of  the  rack 

pitchline  on  the  pitch  circle  of  the  gear  (See  note  below) 

DELAL  =  Act,  the  increment  of  the  involute  generating  angle(See  note 
below) 


CONB 


Initial  value  of  constant  B  of  Lewis  parabola[See  equ.  (E52)] 


DELCON 


AB,  the  increment  of  the  constant  B  of  the 
Lewis  parabola  (See  note  below) 


NRITE  =  Test  constant.  When  NRITE  =  1,  all  the  tooth 
profile  coordinates  stored  in  X2  and  Y2  are 
printed  out.  These  coordinates  are  not  printed 
out  when  NRITE  =  0. 

NOTE:  The  choice  of  ike  magnitudes  of  Aa,  Ai|j  and  AS  is  discussed  in 

Section  3. 

c _ Computations 

Starting  with  line  29  of  the  program  each  location 
of  the  matrix  X  is  filled  with  the  value  -10.00.  This  is  necessary, 
since  the  subsequent  sorting  process  must  determine  the  largest 
value  in  rows  whose  locations  are  either  unfilled  or  contain 
negative  numbers.  If  such  an  unfilled  location  contains  a  zero, 
it  will  be  found  as  the  answer  rather  than  the  actual  x-coordinate 
of  the  left  root  profile. 

Subsequent  to  line  33  various  gear  parameters ,  such  as 
the  pitch  radius  Rp,  the  base  radius  R^  and  the  outside  radius  RQ 
are  computed.  Further,  the  following  quantities  are  determined: 

BEFF 

B  =  b  C(:  =  — nrc~  »  the  actual  distance  to  the  effective 

eff  PD 

point  for  the  given  diametral  pitch 
THETAB  =  0^,  see  equ.  (Al) 
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....... 


The  later  is  computed  with  the  help  of  the 
involutomery  formula  given  by  equ.  (E16)  of  Appendix  E. 
For  this  case: 

T1  Tcs  ’  R1  =  Rp  *  *1  =  6 

and 

T2  =  Tb  '  R2  =  «b  *  4>2  =  0 


so  that 


R  cos  0 

Tb  =  Tcsu -  -  2Rb  [  INV  (0)  ‘  INV  (6)1 


or 


Tb  =  Tcs  cos0  +  2Rb  ^  tan  ~  0  1 


and  finally, 


T  cos8 
cs 


+  2(  tan0 


e  ) 


(2) 


-17- 


Following  line  41  one  finds: 


EPS 

DELTA  = 

BALLOW  = 


b  4.  0  for  use  in  ecju’s.  (B8)  and  (B9) 

IT 

beff  _  tan6,  see  equ.  (B22) 
sin0 

b  ,  according  to  equ.  (Bl) 

effmax 


Lines  45-47  represent  a  test  which  decides  with  the  help  of 
equ.  (Bl)  whether  or  not  the  tooth  is  undercut. 

Lines  48-60  are  applicable  for  undercut  teeth.  The  following 
computations  are  performed: 

A  =  a,  according  to  equ.  (B2) 

TAU  =  EPS,  from  line  42 

The  inner  form  radius  RF  =  Rf  is  obtained  with  the  help  of  equ. 

(D  9) .  This  non-linear  algebraic  equation  is  solved  by  subroutine 

DF.KKER  .  The  radii  R0  and  Rb  represent  the  limit  within  which  the 

root  Rf  is  determined. 

ALPHIN  =  cu n ,  according  to  equ.  (D23) 

ALPHF  =  df in .  according  to  equ.  (D24) 

GAMT  =  y  according  to  equ.  (D16) 


-18- 


PSIIN  =  ^in’  according  to  equ.  (D18) 

PSIF  =  ^fin’  according  to  equ.  (D14) 

Lines  62  to  82  are  applicable  to  teeth  which  are  not  under¬ 
cut.  The  essential  computations  are  given  by: 

A  =  a,  according  to  equ.  (BIO) 

TAU  =  EPS  +  DELTA,  from  lines  42  and  43 [This  ip  n  in  equ's .  (B22),(B24) 

and  (B2b)} 

RF  =  R^,  according  to  equ.  (D12) 

GAMT  =  Ytot’  according  to  e<3u-  (D16)  with  the  above  RF 

ALPHIN  =  ain,  according  to  equ.  (D23) .  The  test  for  a^n 

has  the  purpose  of  avoiding  an  almost  zero  number 
when  a  zero  should  be  obtained 

ALPHF  =  afin>  according  to  equ. (D24)  with  the  RF  of  line  66 

PSIIN  =  0,  according  to  equ.  (D20) 

PSIF  =  in >  according  to  equ.  (D19) 
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The  computations  represented  by  lines  82  to  88  are 
applicable  to  both  types  of  teeth.  In  addition  to  converting 
certain  angles  from  radians  to  degrees  for  printing  out,  the 
following  parameters  of  the  center  of  the  rack  fillet  are 
computed: 

A2  =  a2>  according  to  equ.  (Cl) 

B2  =  b2>  according  to  equ.  (C2) 


Lines  89  to  104  concern  themselves  with  the  determina¬ 
tion  of  the  coordinates  of  the  involute  portion  of  the  tooth.  It 

is  at  first  most  convenient  to  store  XI  =  x.  and  YI  =  y.  ,  which 

inv  7 inv 

are  obtained  from  equations  (A5)  and  (A6) ,  respectively,  in  the 
vectors  XI (I)  and  Yl(l).  For  1=1,  the  angle  a  =  ain-  When  this 
counter  has  attained  its  maximum  value,  the  angle  a  is  just  a 
little  larger  than  af^n-  Note  also  that  this  maximum  value  of  the 
counter  I  is.,  stored  as  IFREV  for  future  use.  For  purposes  of 
subsequent  computations  it  is  desirable  to  have  the  counter  of  the 
coordinate  locations  of  the  tooth  profile  start  at  the  tip  of  the 
tooth  and  end  at  that  root  location  which  corresponds  to 
For  this  reason,  the  sequence  of  the  involute  coordinate  counter 
is  inverted  in  the  operations  of  lines  ICO  to  103.  These  coordinates 
are  now  stored  in  the  vectors  X2  and  Y2.  (Note  that  the  counter 
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number  corresponding  to  IPREV  now  denotes  the  location  of 
the  involute  coordinates  associated  with  the  angle  a^n) . 

The  determination  of  the  coordinates  of  the  tooth 
root,  i.e.  the  trochoid,  starts  in  line  105.  Initially,  the 
vector  Y(I)  and  the  matrix  X(I,J)  are  created  in  the  manner 
described  in  Appendix  C-2.  Subsequently,  a  search  procedure 
determines  the  largest  value  of  the  x-coordinate  X2  associated 
with  a  given  y-coordinate  Y2.  It  is  to  be  recalled  that,  the 
parameters  appropriate  to  the  type  of  tooth  under  consideration 
have  already  been  stored  and  that,  therefore  the  same  basic 
expressions  can  be  used  in  the  program  whether  the  tooth  is  under¬ 
cut  or  not. 


The  following  explains  the  manner  in  which  the  program 
fills  the  vector  Y(I)  and  the  matrix  X(I,J): 

(1)  The  angle  is  made  equal  to  and  the  counter  I, 

for  che  y-coordinate,  and  the  counter  J,  for  the 
number  of  the  circle,  are  initialized. 

(2)  Starting  with  location  25  (line  111),  the  fillet  center 
coordinates  XC  =  x£  and  YC  =  yc  are  computed  according 
to  equ’s.  (C3)  and  (C4)  or  equations  (C5)  and  (C6) . 
Further,  the  y-coordinate  YT  =  y  of  the  effective  noint 
of  the  rack  is  obtained  with  the  help  of  equ.  (B9)  or 
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equ.  (B25).  If  these  computations  were  made  for  the 
first  position  of  the  fillet  circle,  i.e.  for  J  =  1  (in 
the  program  I  =  1  is  used) ,  control  is 
transferred  to  location  80. 

(3)  Starting  with  location  80,  the  value  for  YT  is  multi¬ 
plied  by  1000  and  the  result  is  truncated  such  that  the 
integer  IYT  results.  This  operation,  which  in  essence 
multiplies  a  y- coordinate  given  to  one  one  thousandth  of 
an  inch  by  1000,  is  necessary  for  subsequent  identification 
purposes.  In  addition,  IYT  is  stored  in  vector  Y(I)  in 
location  90.  (Thus,  rather  than  storing  the  actual 
y-coordinate  its  thousand -fold  value  is  used).  If  the 
the  effective  point  of  the  first  circle  is  involved, I  =  1. 
With  the  above  accomplished,  IYT  is  again  made  into  a 
real  number  and  divided  by  1000  (  giving  the  y-coordinate 
to  the  nearest  one  one  thousandth  of  an  inch)  for  further 
work.  It  is  now  labelled  YTTT  [The  y  of  equ.  (C7)  1.  In 
order  to  be  able  to  stop  the  computations  for  points  on 
a  given  fillet  circle, once 

y  <  yc  -  rc  •  (3) 

the  test  of  lines  127  and  128  has  been  devised. 

The  x-coordinate  associated  with  YTTT  is  then  obtained 
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with  the  help  of  equ.  (C9)  of  Appendix  C  and  stored 
in  X(I,J),  where  J  is  the  circle  number  and  I  is  the 
counter  of  the  y- coordinate.  This  sequence,  which  starts 
on  line  124  and  ends  with  line  132  is  continued  until  the 
limit  of  equ.  (3)  is  reached  for  a  given  circle.  In 
each  step  the  counter  I  is  incremented  by  1,  the  integer 
IYT  is  decreased  by  1  and  stored  in  Y(I),  and  finally  the 
associated  x-coordinate  is  stored  in  the  appropriate  X(I,J) 
When  the  computations  for  any  fillet  circle  are  completed, 
control  is  transferred  to  location  150. 

(4)  Starting  with  location  150,  the  stage  is  set  to  compute  the 

x  and  y  coordinates  of  points  on  subsequent  circles :  The 
angle  ip  is  decremented  by  Aip,  and  a  test  is  made  to  check 
whether  has  been  reached.  In  addition,  the  circle 

counter  J  is  incremented  by  1  as  long  as  \p  <  i ^  .  To 
obtain  the  data  for  a  new  circle,  control  is  transferred  to 
location  25  again.  In  case  the  limiting  angle  has 

been  reached,  the  value  J  corresponding  to  the  last  circle 
is  stored  as  JFIN ,  and  control  is  given  to  location  500  to 
commence  the  search  for  the  largest  value  of  X(I,J)  in  a 
given  row  of  the  array. 

(5)  If  control  was  transferred  to  location  25,  the  steps 
described  in  (2) above  are  repeated  for  the  new  position 

of  the  fillet  circle.  In  order  to  be  able  to  form  an  array 
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as  shown  in  Table  C-l  of  Appendix  C,  it  is 
necessary  to  establish  that  location  Y(I)  of  any 
previous  circle  which  contains  an  IYT  of  the  same 
magnitude  as  that  which  corresponds  to  the  effective 
point  of  the  new  circle. 

This  is  accomplished  by  the  sorting  procedure  starting 
in  line  115  and  ending  in  line  121:  The  counter  I,  which 
serves  to  call  on  the  already  existing  Y(I),  is  initialized 
and  YT  is  made  into  the  integer  IYT.  The  desired  location 
Y(I)  is  found  by  comparing  the  content  of  each  of  the 
existing  Y(I)  with  the  value  of  the  above  IYT.  When  the 
test 

Y(I)  -  IYT  <  0  (4) 

is  met,  control  is  transferred  to  location  90  and  the 
present  values  of  I  and  IYT  are  reestablished  as  the 
counter  and  the  content  of  this  Y(I). 

At  this  point  all  operations  for  a  given  circle  may  be 
performed  as  described  beginning  with  the  second  paragraph 
of  section  (3)  above. 

(6)  As  stated  earlier,  once  all  the  possible  Y(I)  and  X(I,J) 
are  filled,  control  is  transferred  to  location  500  to 
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search  for  the  x-coordinates  of  the  root 
profile. 


The  preparation  for  the  determination  of  the  largest 
value  of  X(I,J)  corresponding  to  any  y-coordinate  Y(I)  starts  in 
location  500  (line  139).  YFIN  ,  the  smallest  value  of  the 
y-coordinate  for  which  X(I,J)  has  been  found  earlier,  is 
obtained  by  way  of  YCFIN.  This  variable  is  obtained  with  the 
help  of  the  fillet  center  expressions  (C4)  or  (C6) ,  evaluated  at 
ip  =  •  Then,  similar  to  equ.  (3): 

YFIN  =  YCFIN  -  RC  (5) 

Subsequently,  the  following  operations  are  performed: 

(1) .  The  counter  I  is  initialized.  It  is  used  to  keep 

track  of  the  Y(I)  and  X(I,J)  of  before. 

(2) .  The  counter  K  =  I  4-  I  HlEV  is  introduced  in  location 

501.  This  counter  makes  it  possible  to  add  the 
final  values  of  the  root  coordinates  to  vectors 
X2  and  Y2,  which  presently  only  contain  the  involute 
coordinates  in  locations  from  I  =  1  to  I  =  IFREV. 
(See  line  98). 


(3) .  KFIN  =  K  is  introduced,  in  order  to  be  able 

to  label  the  last  values  of  X2  and  Y2. 

(4)  .  The  sorting  process  for  each  row  of  the  array 

X(I,J)  starts  in  line  144  with  the  initialization 
of  X2(K)  -  -10.00.  The  loop  between  line  145  and 
147  starts  by  comparing  the  above  value  with  that  of 
the  content  of  X(l,l)  i.e.  the  first  value  in  the 
first  row.  If  X(l,l)  is  larger  than  -10.00,  X2(K) 
is  made  equal  to  X(l,l).  If  it  contains  the 
initially  read-in  value  of  -10.00,  X2(K)  remains  as 
-10.00  (See  line  31).  In  the  next  comparison,  the 
content  of  X(l,2)  is  compared  with  the  current  value 
of  X2(K).  In  case  X(l,2)  >  X2(K)  ,  this  value  of 
X(l,2)  is  made  the  current  value  of  X2(K).  This 
sorting  process  is  repeated  until  J  =  JFIN  ,  the 
number  of  the  last  fillet  circle  (See  line  137). 

(5) .  The  outer  loop  of  the  sorting  process  includes 

lines  142  to  151.  Once  the  largest  value  of  any 
row  of  X(I,J)  has  been  stored  in  X2(K),  the  content 
of  the  associated  Y(I)  is  divided  by  1000.00  and 
stored  in  Y2(K).  This  operation  is  repeated  until 
the  last  value  of  Y2  becomes  eoual  to  or  smaller 
than  the  previously  determined  YFIN  . 
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Control  is  transferred  to  location  21,  whenever 
the  vectors  X2  and  Y2  contain  a  complete  set  of  coordinates  for 
the  left  tooth  profile. 

The  determination  of  the  various  parameters,  needed 
for  the  computation  of  the  Lewis  factor,  starts  in  location  21 
with  the  following: 


PHIRO 

-e- 

o 

according 

to  equ.  (E22) 

BETA 

= 

6 

according 

to  equ's.  (F.21) 

THETA 0 

= 

0R 

Ro 

,  see  equ. 

( 029) 

RV 

= 

Rv  • 

according 

to  equ.  ( F28) 

The  search  for  the  point  of  tangency  between  the 
Lewis  parabola  and  the  tooth  profile  is  conducted  according  to 
the  outline  in  Appendix  E-4.  It  is  first  necessary  to  find  the 
value  of  the  counter  K,  used  in  X2(K)  and  Y2(K),  which  corresponds 
to  that  of  the  y-coordinate  RV  of  point  V  of  Figure  E-4.  The 
range  between  this  point  and  the  already  known  YFIN  ,  for  which 
K  =  KFIN  ,  is  the  range  in  which  the  tangency  will  occur. 

For  computational  reasons,  the  point  START  rather 
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than  point  V  is  chosen.  Thus, 

START  =  RV  -  .005  (inches)  (6) 

The  location  of  that  Y2  which  corresponds  to  START 
is  found  in  the  loop  between  lines  165  and  169,  Once  the  number 
of  this  location  is  found,  it  is  labelled  BEGIN  . 

The  determination  of  the  tangency  point  is  carried 
out  in  the  loop  between  lines  172  and  179.  Using  the  counter  J, 
the  value  of  K  is  varied  between 

K  =  BEGIN  +  J  -  1  (7) 

and 

k  =  begin  +  kfin  -  i  ,  (8) 

where  KFIN  is  the  difference  between  the  KFIN  of  line  143  and 
BEGIN  of  line  17  0. 


The  operation  consists  of  the  following  steps: 


(1).  The  value  of  XPAR 


par 


is  computed  according  to 
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equ.  (E32),  using  the  initial  value  of  CCNB  =  B 
[See  equ.  (E31)],  for  the  full  range  of  values  of 
Y2(K).  In  each  step,  the  difference 


DIFF2  =  | X2 (K)  |  -  XPAR  (9) 

is  computed  [See  also  equ.  (' E33)  ]  .  If  this  difference 
is  less  than  or  equal  to  zero,  the  point  of  tangency 
has  been  found  and  control  is  transferred  to  location 
109  (line  180). 

(2).  If  the  point  of  tangency  cannot  be  found  for  a  given 
value  of  CCNB,  this  value  is  decremented  by  DELCCN 
(see  line  178)  and  equ.  (9)  is  used  to  test  a  new 
series  of  XPAR  values.  This  process  is  continued 
until  the  values  of  X2(K)  and  Y2(K),  associated  with 
the  tangency  point,  are  found. 

Finally,  with  X2(K  )  and  Y2(K  )  found,  EX  =  x  is 

C3n  l  3rl 

computed  according  to  equ.  (1-35)  ,  and  the  Lewis  factor  as  well  as 
the  AGMA  geometry  factor  are  determined  with  the  help  of  equ's.  (Fll) 
and  (F.36)  ,  respectively. 
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Output  Parameters 
(All  lengths  are  given  in  inches  and  all  angles  in  degrees). 

The  output  of  the  program  first  states  the  numerical 
values  of  the  input  parameters  N,  BEFF,  H),  THETAD,  CAFTCS,  RC, 
KADD,  DELAL,  DELPSI,  DELCCN  and  CCN B.  Subsequently,  it  is  stated 
whether  or  not  the  tooth  is  undercut.  In  addition,  the  output 
gives  numerical  values  for  the  following  quantities: 


RP 

R 

P 

N 

,  the  pitch  radius  of  the  gear 

2  Pd 

RB 

= 

R  cos  0 

P 

,  the  base  radius  of  the  gear 

RO 

= 

R 

P 

f  KADD 
Pd 

,  the  outside  radius  of  the  gear 

RF 

= 

Rf 

,  the  inner  form  radius,  either  according  to 

equ 

.  (D9)  , 

or  according  to  equ.  (D12) 

EPS 

= 

°b 

2 

+  e  , 

in  degrees 

THETAB 

= 

°b 

,  according  to  equ.  (2) ,  in  degrees 

DELTA 

bef  f 

tan0  ,  according  to  equ.  (B22) ,  in 

Rb 

sine 

degrees 
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TAU 

= 

EPS 

or  EPS  +  DELTA 

GAMT 

= 

^tot 

,  according  to  equ. 

(D16) 

ALFHIN 

= 

a . 
in 

,  according  to  equ. 

(D23) 

ALPHFIN 

= 

. 

fin 

,  according  to  equ. 

(D24) 

PSIIN 

= 

. 

m 

,  according  to  equ. 

(D18)  or  equ.  (D20) 

FSIFIN 

= 

♦fin 

,  according  to  equ. 

(D14)  or  equ.  (D19) 

PHIRO 

= 

*Ro 

,  according  to  equ. 

(H22) 

BETA 

= 

6  , 

according  to  equ’s. 

(E21) 

CAPTO 

= 

To  ’ 

tooth  thickness  at 

Rq ,  computed  according  to 

equ . 

(E 16) 

THETAO 

= 

0R 

0 

,  computed  according 

to  equation  (E15)  and 

furnishing  a  result  identical  to  that  of  equ.  (E29) 
RV  =  R  ,  according  to  equ.  (E28) 
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EX  =  x  ,  according  to  equ.  (E35) 

YLHWIS  =  Y  ,  the  Lewis  factor,  according  to  equ.  (Ell) 

YAGMA  =  ^AGMA  ’  accor<^inS  t0  equ.  (E36) 

TANGENCY  K  =  value  of  counter  K  at  point  of  tangency  between 

Lewis  parabola  and  tooth  profile 
Y2(KTAN)  =  y-coordinate  of  tooth  profile  associated  with 

point  of  tangency 

X2(KTAN)  =  x-coordinate  of  tooth  profile  associated  with 

point  of  tangency 

XPAR  TAN  =  x-coordinate  of  parabola  associated  with  point 

of  tangency 

CONB  =  final  value  of  constant  B  of  Lewis  parabola 

NRITE  =  Test  constant  (see  section  b  above) 

If  NRTTF  =  1,  the  program  causes  the  printing  of  the 
tooth  profile  coordinates  Y2(K)  and  X2(K)  for  all  values  of  K, 
starting  from  the  tip  of  the  tooth  down  towards  its  root.  The 
transition  from  the  involute  to  the  trochoidal  part  of  the  profile 


occurs  where  the  values  of  Y2(K)  are  given  to  one  one- thousandth 
of  an  inch.  Recall  that,  since  the  program  computes  the  left  tooth 
profile,  the  values  of  X2(K)  are  negative. 

II.  Program  LF.WIS  ENVELOPE  (See  Appendix  H) 
a.  Dimensioning  of  Arrays 


Again,  the  vectors  XI  and  Y1  are  used  for  the  initial 
storage  of  the  involute  coordinates  of  the  gear  tooth.  The  sizes 
of  these  vectors  are  also  determined  with  the  help  of  equ.  (1),  or 
by  way  of  program  MATRIX  (see  Appendix  F)  . 

The  vectors  X2  and  Y2  must  again  hold  the  values  of  the 
involute  coordinates  as  well  as  those  of  the  root  coordinates.  The 
required  dimensions  of  these  arrays  may  he  obtained  by  adding  the 
number  of  rows  of  vector  Y1  to  the  following  number  of  trochoid  compu 
tation  intervals.* 

*in  '  ^f in 

Nt  = - -  (see  below)  (10) 

where 

ip-n  =  initial  trochoid  generating  angle  [For  undercut  teeth 

use  equ.  (D18),  for  non-undercut  teeth  equ.  (D20)  is 

*  The  number  of  these  intervals  is  identical  with  the  column  number 
for  the  tooth  root  interval ,  as  given  in  the  output  of  program 
MATRIX  in  Appendix  F.  See  also  equ.  (F8). 
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applicable.  is  positive  or  zero.] 

=  final  trochoid  generating  angle.  [For  undercut  teeth 
use  equ.  (D14)  ,  for  non-undercut  teeth  use  equ.  (D19). 

is  always  negative.  This  explains  the  form  of 
equ.  (10).] 

Avp  =  increment  of  roll  angle  of  the  rack  pitch  line 

(See  earlier.) 

h .  I nput  P a r am e t er s 

The  input  parameters  for  program  LEWIS  ENVELOPE  are  the 
same  as  those  for  program  LEWIS  CIRCLES. 

c  .  Comput at  i ojn s 

With  the  exception  that  MATRIX  X  does  not  exist  in  this 
program,  LEWIS  ENVELOPE  is  parallel  to  LEWIS  CIRCLES  up  to  line 
98:  Various  initial  gear  parameters  are  computed.  It  is  decided 

whether  or  not  the  gear  tooth  is  undercut  and  the  associated  para¬ 
meters  are  determined.  The  involute  coordinates  are  computed  and 
stored  in  vectors  XI  and  Yl,  with  the  last  computation  corresponding 
to  the  outside  radius  of  the  tooth.  In  order  to  he  able  to  obtain 
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tooth  profile  coordinates  which  start  with  the  tip  of  the  tooth, 
the  order  of  the  vectors  XI  and  Y1  is  again  reversed  and  the 
result  is  stored  in  the  vectors  X2  and  Y2.  The  location  number 
of  the  last  pair  of  involute  coordinates  (now  associated  with  the 
angle  oun)  corresponds  again  to  the  value  of  IPREV. 

The  determination  of  the  tooth  root  coordinates  starts 
in  line  99  and  proceeds  in  the  following  manner: 

(1)  The  angle  tp  is  made  equal  to  tp ^ n ,  and  the  counter  I 
is  initialized  at 

I  =  IPRF.V  +1  (11) 

(2)  A  loop,  starting  with  line  104  and  ending  with  line 

112,  computes  the  coordinates  XF.NV  =  x  and  YF.NV  =  v 

env  ■ env 

according  to  equ’s.  (C21)  and  (C22),  respectively,  until 
ip  =  (PSIFD).  The  results  are  stored  in  vectors  X2 

and  Y2 ,  between  the  locations  I  =  IPRF.V  +  1  and  I  =  KFIN. 

The  computations  for  the  LEWIS  and  AGMA  factors  start  in 
line  114  (statement  no.  21)  and  are  identical  to  those  of  program 
LF.WTS  CIRCLES  up  to  and  including  line  151. 

d .  Output  P arameter s 

The  output  parameters  of  program  LEWIS  ENVELOPE  are  identical 
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to  those  of  program  LEWIS  CIRCLES.  The  print-out  of  the  tooth 
profile  coordinates  differs  from  that  in  program  LEWIS  CIRCLES 
somewhat.  For  the  present  case,  it  was  possible  to  separate  the 
involute  coordinates  from  the  trochoid  coordinates  by  appropriate 
headings,  i.e.  INVOLUTE  PLOT  ar.d  TROCHIOD  PLOT. 

3.  SAMPLE  RUNS  OF  PROGRAM  LEWIS  CIRCLES  AND  LEWIS  ENVELOPE 

This  section  first  discusses  the  results  of  two  identi¬ 
cal  sample  runs  which  were  made  with  each  of  the  two  programs.  The 
first  type  of  run  involves  a  non-undercut  36  tooth  gear.  The  second 

concerns  itself  with  an  undercut  12  tooth  pinion. 

Appendix  G  lists  program  LEWIS  CIRCLES  on  pages  G-2  to 
0-9.  The  associated  output  for  the  36  tooth  gear  is  given  on  pages 
G-10  to  G-18,  while  that  for  the  12  tooth  pinion  covers  pages  G-19  to 
G-  37. 

Appendix  H  lists  program  LEWIS  ENVELOPE  on  pages  H-2  to 
H-9.  The  associated  output  for  the  36  tooth  gear  is  given  on  pages 

H-10  to  H-14.  There  are  two  outputs  for  the  12  tooth  pinion.  The  one 

with  the  larger  angle  increment  Aip  is  given  on  pages  H-15  to  H-19,  while 
that  for  the  smaller  increment  A  ip  is  reproduced  on  pages  H-20  to  H-32. 

Finally,  the  results  of  certain  selected  computations  of  LEWIS  and 

AGMA  factors  are  given  in  tabular  form. 
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a.  36  Tooth  Gear 


out  on 

N 

BEFF 

pn 

THETAD 

CAPTCS 

RC 

KADD 

DF.LAL 

DELPSI 

CONB 

DELCON 


1 .  Input  Parameters  [See  also  Section  2-Ib) 

With  one  exception,  all  input  data  of  this  run  are  printed 
top  of  the  first  output  page  (See  p.  G-10). 

=  36  teeth 

=  1.215  in. 

=  1.0,  diametral  pitch 

=  20  degrees 

=  1 . 846  in . 

=  . 04  in. 

=  1.425  in. 

=  .250  degrees 

=  .125  degrees 

=  4.0 

=  .001 
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Further,  on  the  bottom  of  the  first  output  page: 


NRITF  =  1 


2 .  Array  Dimensions 

The  output  of  the  program  MATRIX  on  p.  F-12  gives  the 
following  minimum  dimensions  for  the  various  arrays  associated  with 
the  3b  tooth  gear: 

Length  of  rows  for  involute 

vectors  XI  and  Y1 :  94 

Length  of  rows  and  columns 

for  tooth  root  arrays:  ROW  =  370,  COL  =  85 

This  leads  to  the  minimum  number 

of  locations  for  vectors  X2  and  Y2:  94  +  370  =  464 

As  a  consequence  of  the  above,  the  following  dimension 
statements  are  included  in  the  program  of  this  run  (see  p.  G-2): 

XI  (100)  ,  Y1 (100) 

Y (400)  ,  X (400 ,  100) 

X 2  ( 5 0 0 ) ,  Y2  ( 500) 


-  38- 


1 


3.  Output  Data  (See  also  Section  2-Id) 

The  first  output  statement  on 
fact  that  the  "TOOTH  IS  NOT  UNDERCUT", 
block  is  useful  for  checking  purposes: 


RP 

- 

18.000  in. 

RR 

= 

16.91447  in. 

RO 

= 

19.425  in. 

RF 

17.11373  in. 

THETAR 

= 

7.58392  degrees 

GAMT 

= 

0 . 0  degrees 

EPS 

= 

23. 79196  degrees 

DELTA 

= 

-  8 . 82052  degrees 

TAU 

= 

14.97144  degrees 

AI.PHJN 

= 

8.82052  degrees 

ALPIIFIN 

= 

32.35465  degrees 

PS  1 1 N 

= 

0 . 0  degrees 

PSIFIN 

= 

-10.6258  degrees 

page  G-10  conveys  the 
The  subsequent  data 
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PSIRO 


29.4532  degree?  (This  is  actually  PH1RO  acc.  to 


BETA 

CAPTO 

THF.TAO 

RV 

The 

G  -  1  0  deals  w i 

PX 

Y LEW IS 

YAGMA 

The 

which  arc  of 

TANGF.NCY  K 

Y2  (KTAN) 


Section  2-Id) 


28 . 56269  degrees 
.60383  in. 

.89052  degrees 
19.25831  in. 

first  line  of  the  lowest  block  of  results  on  page 
th  the  essential  output  of  the  program,  i.e.: 

.55278  in.,  the  Lewis  dimension  x 

.3685,  the  Lewis  factor 

.36221,  the  AGMA  geometry  factor 

second  line  of  this  data  block  gives  further  values 
interest  for  checking  purposes: 

=  237,  i.e  the  point  of  tangency  between  the  tooth 

profile  and  the  Lewis  parabola  was  found 
by  the  program  to  occur  at  Y2(237) 

=  1C. 937  in.  ,  y-coordinate  of  tooth  profile 

at  point  of  tangency. 
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X2(KTAN)  =  -1.13277  in.,  x-coordinate  of  tooth  profile 

at  point  tangency 

XPAR  TAN  =  1.132784  in.,  absolute  value  of  x-coordinate 

of  Lewis  parabola  at  point  of 
tangency 

CONB  =  1.809,  final  value  of  constant  B  of  Lewis 

parabola 

The  listing  of  the  coordinates  of  the  left  tooth  profile, 
i.e.  of  the  values  stored  in  the  associated  vectors  X2  and  Y2, 
begins  on  top  of  page  (1-11.  Recall  that  the  values  of  these 
coordinates  are  printed  out  in  such  a  way  that  they  start  at  the 
tip  of  the  tooth  and  end  at  its  root.  The  involute  portion  of  the 
profile  begins  with  K  =  1  and  ends  with  K  =  96.  While  the  first 
set  of  values  of  the  trochoidal  portion  of  the  tooth  profile  at  k  =  97 
is  not  specifically  marked,  it  may  be  recognized  by  the  fact  that 
Y2(97)  represents  the  first  y-coordinate  which  is  given  in  terms  of 
one  one- thousandth  of  an  inch.  This  trochoidal  portion  of  the  profile 
ends  with  K  =  460.  Note  that  the  point  of  tangency  at  K  =  237,  is 
well  above  the  last  point  of  the  profile. 

b.  12  Tooth  Pinion 


(Appendix  G  shows  only  the  output  of  this  run.) 


1. 


With  the  exception  of  the  NRITF  statement,  all  input 
data  are  printed  again  on  top  of  page  G-19: 

N  =  12 

BF.FF  =  1.0526  in. 

PD  =  1 

TUF'TAD  =  20  degree? 

CAFTCS  =  1.57079  in. 

RC  =  . 300  i n . 

KADD  =  1.000  in. 

DF.I.AF  =  .250  degrees 

DFl.PSI  =  .500  degrees 

CONB  =  4.0 

DPI. CON  =  .001 

Further,  as  shown  on  the  botton  of  page  G-19: 

NR  1 TF.  =  1 

2.  Array  D imens ions 

(The  program  listing  containing  the  dimension  statement 
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given  below  is  not  reproduced  in  this  report.) 


The  output  of  program  MATRIX  on  page  F-13  gives  the 
following  minimum  array  dimensions  for  the  given  pinion  data 
(For  rationale,  see  parallel  section  for  the  36  tooth  gear.). 

XI  (150)  ,  Y1(1S0) 

Y (995)  ,  X (995  ,  52) 

X2 ( 1 1 45)  ,  Y2  (1145) 

3 .  Output  Data 

The  first  output  statement  on  p.  G-19  indicates  that 
the  "TOOTH  IF  UNDERCUT".  The  subsequent  data  block  may  be 

interpreted  as  discussed  in  section  2-Id  and  is  parallel  to  the 
one  shown  for  the  36  tooth  gear.  The  first  line  of  the  lowest 
block  of  results  again  represents  the  essential  output  of  the 
run ,  i . e .  : 

F.X  =  .34446  in.,  the  Lewis  dimension  x 

YLF.WIS  =  .  22964  ,  the  Lewis  factor 

YAGMA  =  .2364  ,  the  AGMA  geometry  factor 

The  second  line  of  this  data  block  gives: 


TANGF.NCY  K 


729 


Y2 ( KTAN ) 


5.019  in. 


X2 (KTAN)  =  -  .780016  in. 

XPAR  TAN  =  . 78003  in. 

CONR  =  2.903 


The  listing  of  the  tooth  profile  starts  on  p.  G-20. 

The  involute  portion  of  the  tooth  profile  begins  with  K  =  1  and 
ends  with  K  =  152.  The  trochoidal  part  of  the  profile  extends 
from  K  =  153  to  K  =  1133.  Note  that  the  point  of  tangency,  at 
K  =  729,  is  well  above  the  last  computed  profile  point. 

I  T  .  Sample  Run s_witji  Program  LEWIS  ENVELOPE 
a,  36  Toot h  Gear 

1 .  Input  Parameters 

With  the  exception  of  CONB  =  5.0,  the  input  parameters 
of  this  run  with  program  LEWIS  ENVELOPE  are  identical  with  those 
used  in  the  run  with  program  LEWIS  CIRCLES  (See  section  3-Ia).  As 
before,  these  data  are  listed  principally  on  top  of  the  first  output 
page  of  the  program.  This  output  page  is  reproduced  on  p.  H-10  of 
this  report.  (Note  that  the  format  of  the  first  output  page  of 
program  LEWIS  ENVELOPE  is  the  same  as  that  of  the  first  output  page 


of  program  LEWIS  CIRCLES.) 

2 .  Array  Dimensions 

A  certain  part  of  the  results  of  program  MATRIX  may 
be  used  to  determine  the  dimensions  of  the  vectors  XI  and  Y1  as 
well  as  those  of  the  vectors  X2  and  Y2.  Thus,  from  p.  F - 1 2 : 

Length  of  rows  for  involute  vectors  XI  and  Y1 :  94 

Since  there  is  now  one  set  of  root  profile  coordinates  for  each 
position  of  the  angle  the  number  of  rows  of  these  coordinates 
is  equal  to  the  number  of  columns  given  by  program  MATRIX  under 
the  heading  of  "Number  of  rows  and  columns  for  tooth  root  arrays". 
[See  equ's.  (10)  and(F8)].  Thus,  the  number  of  root  coordinate 
pairs  equals  85  according  to  p.  F-12. 

The  minimum  dimensions  for  X2  and  Y2  then  become  : 

94  +  85  =  179.  Actually,  the  following  larger  dimensions  were 
included  into  the  program  for  this  run  (see  program  LEWIS 
ENVELOPE  on  p.  H-3) : 

XI  (1 80)  ,  Y1 ( 180) 

X2  (1000) ,  Y2  ( 1 000) 
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3.  Output  Data 

The  central  output  data  block,  starting  with  "TOOTH 
IS  NOT  UNDERCUT"  and  ending  with  RV  =  19.25831  is  identical  with 

that  obtained  by  way  of  program  LEWIS  CIRCLES,  as  shown  on 
p.  G-10.  The  first  line  of  the  lowest  output  data  block  again 
furnishes  the  principal  results  of  the  run.  (Recall  that  the 
root  profile  coordinates  are  computed  in  a  different  manner  in 
program  LEWIS  ENVELOPE.)  Thus, 


EX 

= 

.55277  in. 

YLEWIS 

= 

. 36851 

YAGMA 

= 

. 36221 

These  results  are  identical  with  those  obtained  by  way 
of  program  LEWIS  CIRCLES.  The  second  line  of  this  lowest  block 
is  useful  for  comparison  purposes  (see  Section  4): 

TANGENCY  K  =  117 

Y2  (KTAN)  =  16.93642  in. 

X2 (KTAN)  =  -1.132901  in. 

XPAR  TAN  =  1.132924  in. 


CONB  =  1.809 
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The  listing  of  the  tooth  profile  coordinates  starts 


on  top  of  page  H-ll.  As  for  the  comparable  run  with  program 
LEWIS  CIRCLES,  the  involute  portion  of  the  profile  is  given 
from  K  =  1  to  K  =  96.  The  root  part  of  the  profile  extends 
from  K  =  97  to  K  =  182. 

b.  12  Tooth  Pinion 


Program  LEWIS  ENVELOPE  was  run  twice  for  the  12  tooth 
pinion.  The  input  data  of  the  first  run  are  identical  with  those 
used  earlier  in  conjunction  with  program  LEWIS  CIRCLES  (See  section 
3-Ib).  In  the  second  run  the  increment  Aip  was  reduced  to  .05 
degrees  in  order  to  observe  the  influence  of  the  more  closely 
spaced  y-coordinates  of  the  root  profile. 

1 .  Input  Parameters  for  Run  with  =  .5  degrees 

The  input  data  for  this  run,  as  shown  on  p.  H-15, 
are  the  same  as  those  given  in  section  3-Ibl. 

2.  Array  Dimensions  for  Run  with  Aip  =  .  5  degrees 

The  minimum  array  dimensions  for  this  run  may  again 
be  obtained  in  the  manner  discussed  in  section  II -a2  above  with  the 
help  of  the  results  of  program  MATRIX  on  p.  F-13. 


-47- 


Since  the  lengths  of  the  rows  of  involute  vectors  XI  and  Y1 
must  be  at  least  ISO, 

XI (150)  ,  Y1 ( 1 50) 

The  number  of  columns  for  the  root  array  X  of  program  LEWIS 
CIRCLES  was  determined  to  be  52.  Therefore, 

X2  (202)  ,  Y2 ( 2 02) 

3.  Output  Data  for  Run  with  Aip  =  .  5  degrees 

The  central  output  data  block  of  p.  H-15  is  iden- 
tical  with  that  obtained  from  program  LEWIS  CIRCLES  (see  page  G-19). 
The  principal  results,  as  given  in  the  lowest  output  block,  are  as 


fol lows : 

EX 

. 34443 

YLEWIS 

. 22962 

YAGMA 

. 23643 

TANGF.NCY  K  = 

177 

Y2 (KTAN) 


5.01112  in. 


X2 (KTAN) 


-.781767  in. 


XPAR  TAN  =  . 781767  in. 

CONB  =  2.903 

These  results  are  essentially  identical  with  those 
obtained  from  program  LEWIS  CIRCLES.  The  difference  in  the 
y-coordinate  of  the  tangency  point  is  due  to  the  different 
manner  of  computing  the  tooth  root  coordinates.  (See  discussion 
in  section  4  below.)  The  listing  of  all  tooth  profile  coordi¬ 
nates  is  given  on  pp.  H-16  to  H-19  of  Appendix  H. 

4 .  Output  Data  for  Run  with  Aip  =  .05  degrees 

The  input/output  data  of  this  run  with  a  reduced 
angular  increment  are  listed  on  pp.  H-20  to  H-32  of  Appendix  H. 
The  principal  results  are  as  follows: 

EX  =  .34447  in. 

YLF.WIS  =  .22965 

YAGMA  =  .23635 

TANGENCY  K  =  385 
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BEFF 
(in.  ) 

CAPTCS 
‘  (in. ) 

RC 

(in .  ) 

KADD 
(in. ) 

1.000 

1 . 57079 

.2  39 

1.000 

\  1.0526  7.57079 

I  it  n 

I 


300  l.oon 


Lewis  Factor 


.  24189 


.25493 


.  33822 


.  36134 


. 39494 


.40209 


.43683 


.45828 


. 3351 7 


. 34831 


. 37016 


. 47858 


. 53370 
_________ 

.  24321 


AGMA 

Factor 


. 24807 


.25922 


. 32876 


. 34784 


. 37542 


. 38121 


.40976 


.42337 


. 35521 


.44417 


.48946 
72364  0 
. 24817 


all  runs:  Diametral  pitch  is  unity  and  hob  pressure  angle  is 
All  runs  made  with  program  LEWIS  CIRCLES.) 


Table  1 

Selected  Computations  of  Lewis  and 
AGMA  Geometry  Factors 
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rr 

i  I 


lL 


Y2 (KTAN) 


X2 (KTAN) 


XPAR  TAN 


CONB 


5 .02859  in. 


-.777906  in. 


777909  in. 


2.903 


III.  Selected  Computations  of  the  Lewis  and  the  AGMA  Geometry 


Factors 


Table  1  shows  the  results  of  selected  determinations 
of  the  Lewis  and  the  AGMA  geometry  factors.  In  all  cases,  the 
computations  were  made  with  the  help  of  program  LEWIS  CIRCLES 
and  a  diametral  pitch  of  one,  together  with  a  20  degree  hob 
pressure  angle  were  used. 


4.  DISCUSSION 


This  section  first  presents  a  discussion  of  certain 
differences  in  the  results  for  the  two  sample  gears  of  section  3, 
as  obtained  with  programs  LEWIS  CIRCLES  and  LEWIS  ENVELOPE.  (Table 
2  juxtaposes  the  essentials  of  these  results.) 


*  m  *m 


Subsequently,  a  number  of  special  features  and 
peculiarities  of  both  programs  are  pointed  out. 

I .  Comparison  of  Results  for  36  Tooth  Gear  and  12  Tooth  Pinion 

When  comparing  the  results  of  the  two  programs  for  a 
given  set  of  data,  it  must  be  kept  in  mind  that  the  point  of 
tangency  of  the  Lewis  parabola  is  generally  located  in  the  root 
portion  of  the  tooth  profile  and  that  the  programs  produce  the 
root  profile  in  different  ways. 


Program 

YLEWIS 

YAGMA 

Y2 (KTAN) 
(inches) 

X2 (KTAN) 

( inches) 

- 

XPAR  TAN 
(inches) 

CONB 

CIRCLES 

Aip  =  .125 

36 

. 36852 

.3622 

16.937 

1.13277 

1 .13278 

1 . 809 

ENVELOPE 

Alp  =  .125 

36 

. 36851 

.36221 

16.9364 

1.13290 

1.13292 

1.809 

CIRCLES 

Aip  =  .  500 

1  2 

. 22964 

.  23640 

5.0190 

0.78002 

0.78003 

2.903 

ENVELOPE 

Aip  =  .500 

12 

i 

.22962 

.23643 

5.0111 

0.78173 

0.78177 

2.903 

ENVELOPE 

Aip  =  .050 

12 

. 22965 

. 23635 

5.0286 

0. 77791 

0. 77791 

2.903 

Table  2 


Comparisons  of  Results  of  Programs 
LEWIS  CIRCLES  and  LEWIS  F.NVELOK. 


Program  LEWIS  CIRCLES  furnishes  one  x-coordinate 
of  the  tooth  root  for  every  one  one-thousandth  of  an  inch  of 


the  y- coord inate .  The  accuracy  of  these  x -coordinates  improves 
as  the  increment  Aip  is  made  smaller,  i.e.  as  more  circles  become 
involved  in  this  simulation  of  the  drafting  process.  (See  Appendix 
C-2.  ) 


Program  LEWIS  ENVELOPE  produces  accurate  x -coordinates 
of  the  root  profile,  because  of  its  analytical  origin.  To  obtain 
these  x-coordinates  in  one  one-thousandth  of  an  inch  increments, 
the  increment  Aip  has  to  be  quite  small. 

Let  the  results  in  Table  2  for  the  3b  tooth  gear  be 
compared  first:  Roth  sets  of  results  were  obtained  with  =  .125, 
and  they  are  essentially  identical.  The  spacing  of  the  v  -  coord i na tes 
due  to  program  LEWIS  ENVELOPE  may  be  seen  on  p.  H-13  of  Appendix  H. 
The  point  of  tangency  is  at  Y2(117)  =  lb. 9364  in.  The  preceding  as 
well  as  the  succeding  /alues  of  the  y- coord inat es  are  approximately 
.006  inches  away.  Because  of  the  magnitude  of  these  distances,  it  is 
probably  fortuitous  that  such  close  agreement  could  be  obtained. 

While  the  results  for  the  Lewis  and  AGMA  factors  as 
well  as  for  the  constant  B  are  also  essentially  identical  in  the 
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case  of  the  12  tooth  pinion,  there  are  some  small  differences 
in  the  x  and  y-coordinates  of  the  root  profiles.  Even  though 
the  increment  AiJj  =  .500  is  quite  large,  it  had  to  be  used  in 
running  program  LEWIS  CIRCLES,  since  the  associated  array  X 
already  had  a  995  x  52  dimension  and  with  that  required  a  very 
large  amount  of  computer  core.  (See  output  of  program  MATRIX 
on  p.  F-13.)  When  this  increment  was  used  to  run  program 
LEWIS  ENVELOPE,  the  magnitude  of  the  spacing  of  the  y-coordinates 
near  the  point  of  tangency  was  approximately  .020  inches.  (See 
K  =  177  on  p.  H - 1 9 . ) 

When  program  LEWIS  ENVELOPE  was  run  with  A</>  =  .05. 
the  spacing  of  the  y- coor idnates  surrounding  the  point  of  tang¬ 
ency  was  of  the  order  of  .002  inches  (See  K  =  385  on  p.  H-27.). 
This  result  appears  to  be  the  best  of  the  three. 

Aside  form  these  small  discrepancies,  both  programs 
furnish  excellent  results. 

1 1  Special  Features  of  Program  LEWIS  CIRCLES 

a.  All  input  data  must  be  expressed  in  terms  of  unity 


-54- 


diametral  pitch. 


The  program  cannot  be  run  for  a  hob  fillet  radius 
RC  =  0,  since  its  simulation  of  the  drafting  process 
requires  such  a  radius.  This  limitation  can  be  over¬ 
come  by  choosing  a  very  small  radius  for  the  case  of 
a  sharp  cornered  hob.  (By  letting  RC  =  .040  in.  for 
the  unity  diametral  pitch  computation,  an  actual  gear 
of  60  diametral  pitch  will  be  cut  by  a  cutter  with  a 
radius  of  .00067  in.) 

The  angular  increment  A  a.  must  be  chosen  sufficiently 
small  if  one  desires  a  fairly  close  spacing  of  the 
y- coordinates  of  the  involute  part  of  the  tooth  profile 
for  drawing  purposes.  A  Aa  =  .25  degrees  caused  a  spac¬ 
ing  of  between  .015  and  .030  inches  for  the  36  tooth 
gear  and  a  spacing  of  between  .003  and  .015  inches  for 
the  12  tooth  pinion.  The  larger  spacing  will  always 
occur  near  the  tip  of  the  tooth. 

The  choice  of  the  root  angular  increment  is  especially 
important.  The  smaller  Ai(i,  the  larger  the  number  of 
circles  for  a  given  y-coordinate  of  the  root  region. 


and  with  that  the  greater  the  accuracy  of  the 
resulting  root  profile.  This  becomes  crucial 
when  the  hob  fillet  radius  RC  is  small.  For 
RC  =  .040  in.  an  increment  =  .125  degrees 
gave  good  results.  With  RC  =  .300  in.,  the 
increment  Ai|i  =  .500  deg.  was  satisfactory. 

(There  is  a  trade-off  between  accuracy  and  core 
space  requirement.) 

Because  the  number  of  overlapping  circles 
decreases  near  the  bottom  of  the  root  profile, 
it  is  a  good  idea  to  check  whether  the  v-coor- 
dinate  of  the  point  of  tangency  is  at  least  one 

quarter  of  the  magnitude  of  the  hob  fillet  radius 
away  from  the  final  y-coordinate  of  the  root 
profile.  If  this  is  not  the  case,  it  is  best 
to  diminish  the  angle  PSIFIN,  i.e.  continue 

further  with  the  drafting  simulation.  (The  above 
serves  only  as  a  cautionary  note,  since  during 

all  trial  runs  the  point  of  tangency  was  well 

away  from  the  end  of  the  tooth  root.) 

Since  the  root  profile  coordinates  are  stored  at 
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.001  inch  intervals  of  the  y-coordinate ,  the 
point  of  tangency  will  be  accurate  as  long  as 
the  root  profile  is  accurate. 


III.  Special  Feature  of  Program  LEWIS  ENVELOPE 

a.  All  input  data  must  also  be  expressed  in  terms 
of  unity  diametral  pitch 

b.  The  program  can  be  run  for  all  values  of  the  hob 
fillet  radius,  including  RC  =  0. 

c.  The  above  remarks  concerning  the  angular  increment 
Aa  also  apply  here. 

d.  While  the  y- coordinates  of  the  root  portion  of  the 
tooth  profile  are  given  at  intervals  of  .001  in. 
by  program  LEWIS  CIRCLES,  the  spacing  of  these 
coordinates  by  program  LEWIS  ENVELOPE  depends 
entirely  on  the  magnitude  of  the  angular  increment 

Aip.  If  this  increment  is  too  large,  the  y-coordinates 
of  the  root  profile  may  have  a  fairly  wide  spacing. 
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Since  the  point  of  tangency  of  the  Lewis 
parabola  can  only  be  found  for  existing 
y- coordinates ,  this  point  will  only  then  be 
accurately  determined  if  the  spacing  of 
these  coordinates  is  of  the  order  of  .001 
inches.  Since  there  is  no  core  problem  with 
program  LEWIS  ENVELOPE,  <_  .05  degrees  can 
read i 1 v  be  used . 
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APPENDIX  A 


DETERMINATION  OF  INVOLUTE  COORDINATES 

Figure  A-l  shows  an  involute  gear  tooth,  originating 
at  its  base  circle,  in  the  standard  vertical  position  as 
defined  for  this  report. 

It  is  assumed  that  the  circular  tooth  thickness  t^  =  W' 
at  the  base  circle  and  the  associated  angle  65  are  known.  To 
determine  the  x  and  y  coordinates  of  the  arbitrary  point  T  on 
the  involute,  the  unit  vectors  nou  and  nut ,  in  directions  OU 
and  UT,  respectively,  must  first  be  determined.  Thus, 

eh  e, 

n  =  cos(90  +  ij—  -  a)  i  +  sin(90  +  *—  -  a)  j  (Al) 

OU  ^  ^ 

where 

a  =  the  roll  angle  of  the  generating  tangent  UT 

t. 

0^  =  g—  ,  where  R^  represents  the  base  radius 

b 

Equation  (Al)  may  be  simplified  to: 

°b  °b  -  (A2) 

nQu  =  -  sin (^-  -  a)  i  +  cos(^  -  a)  j 

The  unit  vector  n  is  found  by  cross  multiplication  of  nQU 
by  the  unit  vector  k: 


.i-l 
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T(xinv’y 


V— < 


Point  V  is 
origin  of  involute 

H 

at  base  circle  _ 

2 


Figure  A-l 

Involute  Coordinates 


(A3) 


w 


n 


ut 


=  E 


x  n 


ou 


=  -  cos (^ —  -  ot)  i  -  sin (75—  -  a)  J 


T~ 


The  vector  5T  may  now  be  written  in  the  following  form: 


(5T  = 


R,  n  +  R,  an  . 
b  ou  D  ut 


(AA) 


Substitution  of  equations  (A2)  and  (A3)  into  the  above  and 
subsequent  separation  into  x  and  y-components ,  gives  the 
following  coordinates  of  an  arbitrary  point  T  on  the  left 
involute  profile: 


6b  0b 

xinV  =  Rb  [  -  sin  (^—  -  a)  -  a  cos(^-  -  a)] 


and 


eb  9b 

yinv  =  Rb  [  cos  (tj-—  -  a)  -  a  sin(^  -  a)] 


(A5) 


(A6) 
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APPENDIX  B 


DETERMINATION  OF  COORDINATES  OF  PATH  OF  EFFECTIVE  RACK  POINT 


In  order  to  be  able  to  plot  the  root  of  a  gear 
tooth,  as  cut  by  a  rack  cutter,  it  is  first  necessary  to 
obtain  the  x  and  y  coordinates  of  the  trochoidal  path  described 
by  the  effective  point  of  the  rack  cutter  tooth. 

Figure  B- 1  shows  the  effective  points  of  rack  cutter  teeth 
with  and  without  fillet  radii  as  well  as  the  associated  effective 


a.  EFFECTIVE  POINT 

FOR  SHARP  CORNERED 
RACK 


b.  EFFECTIVE  POT NT 
OF  RACK  WITH 
FI  FEET  RADIES  r 


Figure  B-  1 

Effective  Point  and  Effective  Distance  bcff  of 
Racks  with  and  without  Fillet  Radii 
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dis*-  trices  be£f.  When  there  is  no  fillet  radius ,  (see 
Fig.  B-la)  the  corner  becomes  the  effective  rack  point 
and  bg£f  is  represented  by  the  distance  from  the  acting 
pitch  line  of  the  rack  to  this  corner.  If  the  rack  has  a 
fillet,  the  point  of  tangency  of  the  rack  flank  and  the 
fillet  circle  becomes  the  effective  point  and  be^^  is 
measured  to  this  point. 

The  magnitude  of  bg^^  decides  for  a  given  gear  or 
pinion  whether  the  teeth  are  undercut.  Figure  B-2 
illustrates  the  applicable  criterion. 

If  the  effective  point  of  the  rack  cutter  tooth 
moves  horizontally  along  a  line  which  passes  through  point 
U(the  origin  of  the  involute)  or  above  it,  the  involute 
tooth  will  not  be  undercut.  Thus  the  limiting  distance 

Vffmax  is  ?Jven  by 

bef fmax  =  Rp  -  Rfe  cosG 

bef fmax  =  RP(1  ‘  C°s20)  =  *P  sin20  (B1) 

Thus,  generally, 

If  begf  <_  RpSin20,  the  gear  tooth  will  not  be  undercut.  The 
effective  point  of  the  rack  initially  contacts  the  gear  tooth- 
to-be-cut  above  its  base  circle. 

If  >  RpSin20,  the  gear  tooth  will  be  undercut.  The 

effective  point  of  the  rack  initially  contacts  the  tooth - 
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Figure  B-2 

Maximum  Allowable  beff  which 


avoids  Undercutting 


! 


to-be-cut  below  its  base  circle. 

The  following  first  gives  derivations  for  the  trochoid 
coordinates  of  the  path  of  the  effective  rack  point  for  both  of  the 
above  cases.  Appendix  ('  shows  how  to  obtain  the  associated 
coordinates  of  the  path  of  the  center  of  the  fillet  circle. 

1_.  _  C0_0RD  1 NAT 1 :  s  PATH  OF  EFFECT  I VI:  RAC KP01  NT 

FOR  b  r f  >  R  s  i  n  ^  0 
e  t  r  p 

Figure  B-3  shows  a  rack  cutter,  which  is  about  to  cut  the 
root  of  the  indicated  tooth,  in  its  standard  position.  The 
flank  of  the  rack  tooth  is  tangent  to  the  initial  involute 
point  V. 

The  angle  tp ,  which  defines  the  rolling  of  the  rack 
pitch  line  on  the  gear  pitch  circle,  is  zero  when  line  OP 
is  at  the  rack  pressure  angle  6  to  the  initial  involute 
tangent  OV. 

As  the  angle  tp  increases  in  a  CCW  direction,  a  part 
of  the  undercut  root  of  the  gear  tooth  is  generated.  A  clock¬ 
wise  decrease  in  the  angle  tp  causes  the  completion  of  this 
root.  In  actual  gear  hobbing,  it  would  also  cause  the 
generation  of  the  involute  portion  of  the  tooth.  The  deter¬ 
mination  of  the  required  maximum  positive  and  negative  values 
of  the  angle  tp,  which  are  necessary  to  form  a  sufficient  por¬ 
tion  of  the  root  profile,  together  with  the  determination  of 
the  inner  form  radius  of  undercut  teeth,  is  shown  in  APPENDIX  C. 
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When  t>e££  >  RpSin20,  the  distance  from  the 
instant  center  P  to  point  R.  becomes  for  the  indicated 
standard  position,  i.e.  for  =  0 

a  =  (  Rp  -  be£^)  tan  B  (B2) 

For  non- zero  values  of  ip .  the  distance  PR  is  given  bv 

PR  =  a  +  Rpi|*  m) 

To  define  the  position  vector  OS,  from  the  center 
of  the  gear  to  the  effective  rack  point  S,  it  is  necessary 
to  define  unit  vectors  for  the  directions  of  lines  OP  and 

_  I- 

PR,  respectively. 

0  6 

nQp  =  cos  (  90  +  ^  +  e  +  ip)  i  +  sin(  q0  +  -^  +  B  +  if' )  J  (B4) 

or 


n  =  -sin(  +  0  +  tfi)  i  +  cos(  —  +  B  +  ij;)  i  (B5) 

op  2  2 

Further  , 

Bh  B, 

npr  =  -k  x  nQp  =  cos  <2—  +  0  +  i|>)  i  +  sin(j-  +  B  +  <|0  j  (B6) 

The  position  vector  OS  may  now  be  defined  as: 

OS  =  Rp  nQp+  (  a  +  Rp’lOfipr  “  bgff  nop 

Substitution  of  equ's  (B5)  and  (36)  and  separation  of  x  and  y 
components  gives  the  following  coordinates  for  the  path  of  the 
effective  rack  point: 
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J 


1 


xtr  =  -  (Rp  -  beff)  sin  +e  +  i^) 

e. 

+  (Rp^  +  a)  cos  - 1-  6  +  ip)  (B8) 

and 

0, 

ytr  =  (Rp  -  beff)  cos  (j2  +  e  +  ip) 

°h 

+  (R  ip  +  a)  sin  +  0  +  ip)  (B9) 

2.  COORDINATES  OF  PATH  OF  EFFECTIVE  RACK  POINT 

FOR  b  ce  <  R  sin20 

_ eff  -  p _ 

Figure  B-4  shows  the  effective  point  of  the  rack  making 
initial  contact  with  point  Q  on  the  (future)  involute  profile. 
(This  is  considered  standard  position  for  the  non-undercut  teeth.) 
This  point  must  be  located  on  the  line  of  action  UP.  Note  that, 
point  V  again  represents  the  initial  point  of  the  involute.  The 
angle  ip ,  which  again  defines  the  rolling  of  the  rack  pitch  line 
on  the  pitch  circle  of  the  gear  is  zero  in  the  indicated  standard 
position,  when  line  OP  makes  the  rack  pressure  angle  0  with  the 
instantaneous  base  radius  OU  of  length  R^ .  A  clockwise  decrease 
in  ip  starts  the  generation  of  both  the  involute  and  root  profiles. 
For  a  discussion  of  the  limits  of  angle  ip  see  APPENDIX  C.  (Let  the 
rack  in  Figure  B-4  be  considered  to  be  in  standard  position.)  For 
the  present  case,  with  beff  _<  Rpsin?0,  the  distance  a  becomes: 
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a  = 


eff 


tan  0 


(BIO) 


This  represents  the  distance  PR  when  =  0 .  In  general 


PR  =  a  +  Rp^ 


(Bll) 


To  define  the  position  vector  5?,  from  the  gear 
center  0  to  the  effective  rack  point  S,  it  is  necessary  to 
introduce  unit  vectors  for  the  directions  of  lines  "DP  and 
PR,  respectively,  Now,  (see  Figure  B-4) : 

h  =  cos(  90  +  n  +  \j>  )i  +  sin(  90  +  n  +  ^  )j  (B12) 

or 


hop  =  -sin(n  +  <l>)i  +  cos(  n  +  ^  )j 
further , 

n  =  -k  x  n  =  cos(  n  +  ^  )i  +  sin(  n  +  ^  )j 
pr  op  J 

The  angle  n  is  defined  as  follows: 

°b  +  . 

n  - 

where 


(B13) 


(B14) 


(B15) 


C  =  6  -  *  VOU  ,  (B16) 

and  G  is  the  pressure  angle  of  the  rack  . 
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Further  , 


*  VOU  =  =  §5  (B17) 

b  b 

by  the  properties  of  the  involute  curve. 

But 

QU  =  UP  -  PR  ,  (B18) 


where 


UP  =  tan0 


(B19) 


and 


PR 


eff 

sin6 


(B20) 


Finally  , 


*  VOU  =  tan0  - 


eff 

R^sinG 


and,  according  to  equ's  (B15)  and  (B16) : 


n  = 


T 


+  o  + 


eff 


-  ta-'G 


RusinG 


The  position  vector  OS  may  now  be  written: 


(B21) 


(B22) 


U5=R  h  +(a  +  RiJ»)n  -  b  ccn 

p  op  pr  pr  eff  op 


(B23) 


Substitution  of  equ's  (B13)  and  (B14)  into  the  above, 
followed  by  the  separation  of  x  and  y  components  leads  to  the 
following  coordinates  of  the  effective  point  for  the  present 


casa 


xtr  ~  "  (  Rp  "  t>eff)sin(n  +  ♦)  +  (Rpt|*  +  a)cos(n  +  4>) 

and 

ytr  =  (  Rp  ‘  beff)cos(n  +  *)  + 

(Rp  +  a)sin(n  +  <P) 


(B24) 


(B25) 
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APPENDIX  C 


GENERATION  OF  TOOTH  ROOT  IN  PRESENCE  OF  FILLET  RADIUS  ON 

RACK  CUTTER 


The  present  appendix  shows  two  methods  of  obtaining 
the  coordinates  of  a  tooth  root,  i.e.  the  portion  of  the  tooth 
below  the  involute  profile,  when  the  rack  cutter  has  a  fillet 
radius  r  .  Now,  the  root  profile  is  not  generated  by  the  effective 
point,  as  is  the  case  for  a  sharp  cornered  rack,  but  it  becomes 
the  envelope  of  the  various  positions  of  the  fillet  circle.  This 
fillet  circle  serves  as  a  generating  curve. 

The  first  method  simulates  the  applicable  drafting 
method  of  finding  an  envelope  by  a  purely  numerical  computer 
procedure.  The  center  of  the  fillet  circle  is  first  found  for 
a  given  position  of  the  rack  pitch  line.  Then,  the  x-coordinates 
of  'Oi nts  on  the  appropriate  portion  of  the  associated  circle 
are  'round  as  functions  of  equally  spaced  y-coordinates .  This 
process  is  repeated  for  other  positions  of  the  rack  pitch  line. 

For  all  angles  ip  the  x-coordinates  of  circle  points,  which  are 
associated  with  certain  y-coordinates  ,  are  stored  in  the  same 
row  of  a  matrix.  If  the  angular  increment  of  pitch  line  rotation 
is  sufficiently  small,  that  x-coordinate  which  is  closest  to 
the  tooth  center  line  represents  the  desired  point  on  the  root 
profile . 
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Fipure  C-l 

Center  of  Fillet  Radius  in 
Standard  Position 


The  second  method  was  first  shown  by  R.  G. 


Mitchiner  and  H.H.  Mabie  in  "The  Determination  of  the  Lev.  s 
Form  Factor  and  the  AGMA  Geometry  Factor  J  for  External  Spur  Gear 
Teeth".  Paper  No.  80-DET-59,  Design  Engineering  Technical 
Conference,  Beverly  Hills,  CA. ,  Sept.  28  -  Oct.  1,  1980.  It 
makes  use  of  the  fact  that,  the  instantaneous  centers  of  curva¬ 
ture  of  a  generating  curve  and  its  envelope  lie  on  a  line 
through  the  instant  center  of  rotation  of  the  two  planes 
involved.  The  contact  point  between  the  generating  curve  and 
the  envelope  also  lies  on  this  line  and  it  represents  a  point  on 
the  final  envelope  profile. 

Before  either  of  these  procedures  can  be  shown  in 
detail,  it  is  necessary  to  find  a  method  for  determining  the 
coordinates  of  the  center  of  the  fillet  radius  as  a  function  of 
the  angle  iK  This  is  accomplished  by  the  appropriate  adaptations 
of  the  expressions  given  in  Appendix  B. 

1.  Coordinates  of  Center  Of  Fillet  Radius  for  Undercut  and 

Non-Undercut  Teeth 


Figure  C-l  shows  a  rack  cutter  tooth  as  it  appears 
both  for  undercut  and  non-undercut  gear  teeth  in  the  standard 
position  of  Figure  B-3  and  B-4,  respectively.  The  distance 
from  point  R,  on  the  rack  pitch  line,  to  the  effective  rack 
point  is  represented  by  The  dist«nce  PR,  from  the  instant 

center  P  to  point  R,  which  corresponds  to  the  effective  rack  point  , 
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is  given  by  length  a.  [  As  before,  either  equ's.  (B2)  or 
(BIO)  are  applicable.]  If  one  wishes  to  express  the  position 


of  the  center  C  of  the  right  hand  fillet  (of  radius  r  )  with 
respect  to  the  instant  center  P  in  this  standard  position, 
one  obtains: 


a2  =  ^2  =  a  “  rccos6  (Cl) 

and 


b2  =  ^  =  bef f  “  rcsinQ 


(C2) 


With  the  above  it  is  now  possible  to  adapt  the 
coordinate  expressions  for  the  effective  rack  point  to  new 
ones  for  the  fillet  center  C  by  replacing  the  distances  a  and 
be^  by  a.^  and  b^  ,  respectively.  Thus,  one  obtains  for  undercut 
teeth,  according  to  equations  (B8)  and  (B9) : 


=  -  (R  -  b9)  sin  (  -z~  +  0  + 


c  "p  J2}  '  2 


(C3) 


+  (Rp'i*  +  +  9  +  4») 

♦ 


C-4 

-78- 


and 


yc  ■  (Rp  -  b2)  cos  (  ~  +  e  +  4,) 

eh 

+  (Rpip  +  a^  sin  (  2~  +  6  +  '(') 

Similarly,  the  coordinates  of  the  fillet  center  C  for 
teeth  become  with  the  help  of  equ's.  (B24)  and  (B25) : 

x£  =  -(R  -  b2)  sin  (n  +  <j0  +  (Rp<Jj  +  a2)  cos  (n  +  <P) 

and 


yc  =  (Rp  -  b2)  cos  (n  +  <10  +  (Rp<l<  +  a2)  sin  (n  +  vp) 
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(C4) 


non-undercut 


(C5) 


(C6) 


2. 


Iterative  Method  of  Obtaining  Root  Profile 


Figure  C-2  indicates  how  successive  fillet  circles 
may  be  used  to  obtain  the  root  profile  envelope  for  undercut  as 
well  as  non-undercut  teeth.  Point  represents  the  fillet 
center  when  the  rack  cutter  is  in  the  usual  standard  position. 
The  operative  pitch  radius  is  OP  =  and  the  coordinates  of 

point  Cp  with  respect  to  point  P,  are  given  by  a ^  and  \>2  (as 
defined  in  the  previous  section).  Centers  C2,  and  result 
after  successive  clockwise  rotations  of  the  rack  pitch  line,  wi 


associated  pitch  radii  R^.  Rp3  and 


respectively . 


Figure  C-3  represents  an  enlarged  view  of  the 
relevant  segments  of  two  such  fillet  radii.  It  will  serve  to 
illustrate  the  following  explanation  of  the  computational 
procedure  which  is  used  to  simulate  the  drafting  method  of 
obtaining  the  root  envelope: 


a.  The  coordinates  xc  and  yc  of  centerpoint  C-^ 

of  the  first  circle  are  computed  with  the 
help  of  cither  eauations  (C3)  and  (CA)  or 
equations  (C5)  and  (C6) ,  depending  on  whether 
the  tooth  is  undercut  or  not. 

The  counter  J  is  used  to  indicate  the  number 
of  the  circle.  J  =  1  for  the  present  case 
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k 


Root  Envelope 


Center  1  (Ci  ) 


■  Center  2  (C.i ) 


Figure  C-2 

Generation  of  Root  Envelope 
by  Fillet  Circles 


Center  3  (C3) 


Center  4  (C4) 


\'  /J 


Effective  Rack  Points  of  Circles 


(Circle  1) 


-  o2  V  X(I,1)  X(I,2) 

(Circle  2)  \\  | 


Fillet  Circle 
has  Radius  r 


Center  line  of  tooth 


-Lowest  Points 
of  Circles 
.  (y  =  yc  -  rc) 


Figure  C-3 

Iterative  Determination  of 
Points  on  Root  Envelope 


and  generally,  for  that  circle  which 
corresponds  to  the  intial  angle  i|k  . 

The  total  number  of  circles  depends  on 
the  increment  Ai|>  as  well  as  on  the  final 
angle  .  (These  initial  and  final 

angles  are  defined  in  Appendix  D) . 

b.  Since  the  center  of  the  fillet  circle 
never  lies  at  the  origin,  the  equation 
for  points  on  the  circle  is  given  by: 

(x  -  xc)2  +  (y  -  yc)2  =  r2  (C7) 

The  x-coordinates  of  points  on  this 
circle,  for  a  given  value  of  the  y-coor- 
dinate,  may  be  obtained  from 

x  =  >  c  +  /  x\  -  (y  -  yc) 2  (C8) 

For  practical  purposes  it  is  only  necessary 
to  obtain  the  x-coordinates  of  points  whose 
y- coordinates  are  located  between  the  effective 
rack  point  of  a  given  circle  and  y  =  yc  -  rc 
(See  Figure  C-3).^*^  Since  the  center  of  all 

'  'Points  above  the  effective  point  on  the  rack  flank  are  not 
involved  in  the  generation  of  the  final  root  profile. 
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generating  circles  lie  in  the  second 
quadrant,  their  x- coordinates  will  always 
be  negative. 

Further,  since  the  circle  segments  of 

interest  will  always  be  located  to  the 
right  of  their  center  points,  the  values 

of  all  x-coordinates  will  be  larger  than 

those  of  their  associated  x  's.  For  this 

c 

reason,  the  positive  sign  must  be  used  in 
equation  (C8) ,  i.e.: 

x  =  xc  +  /  r’  -  (y  -  yc)2  (C9) 

c.  The  computation  is  continued  by  determining 
the  y-coordinate  of  the  effective  rack 
point  of  circle  1  either  with  the  help  of 
equ.  (B9)  or  equ.  (B25) ,  depending  on  whether 
the  tooth  is  undercut  or  not. 

The  y-coordinates  are  assigned  the  counter  I 
and  they  will  be  stored  in  the  vector  Y(I). 

For  this  specific  computation  1=1.  The 
associated  x-coordinate  is  then  found  with 
the  help  of  equ.  (C9) ,  and  its  values  is  stored 
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in  the  matrix  X(I,J).  Since  this  is  the 
first  x-coordinate  of  circle  1,  it  will 
be  located  by  X(l,l).  Subsequently,  the 
first  y-coordinate  is  decremented  by  .001 
inch  and  its  associated  x-coordinate  is 
again  obtained  with  the  help  of  equ.  (C9) . 

The  value  of  the  y-coordinate  is  stored  in 
Y(2),  while  the  x-coordinate  belongs  to 
X(2,l).  This  process  is  repeated  until 
y  =  yc  -  rc  (with  the  yc  of  circle  1) ,  filling 
the  appropriate  Y(I)  and  X(I,1). 

d.  The  angle  is  now  decremented  by  Aip ,  as 
indicated  by  Figure  C-3.  J  =  2  for  this 
second  position  of  the  fillet  circle.  Again, 
the  coordinates  xc  and  yc  of  point  C2  are 
found,  together  with  the  y-coordinate  of  the 
effective  rack  point  of  circle  2.  The  later 
is  identical  to  one  of  the  y-coordinates  of 
circle  1.  In  order  to  find  the  value  of  I 
which  corresponds  to  this  Y(I)  of  circle  1, 
all  previously  computed  values  of  Y(I)  are 
compared  to  the  present  one.  (This  sorting 
process  is  made  possible  by  multiplying  the 
values  of  all  the  y-coordinates  stored  in  Y(I) 
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by  one  thousand  and  then  truncating  them 
to  make  a  four  digit  integer  out  of  them. 

(See  Program  description  in  Section  2  ). 

Once  the  specific  value  of  I  =  f  has  been 
found,  it  is  assigned  to  X(f,2)  for  storing 
the  x-coordinate  of  the  effective  point  of 
circle  2.  The  process  of  decrementing  the 
y- coordinates  of  the  second  circle  by  .001 
inch  each  time  and  computing  the  associated 
x-coordinates ,  until  the  lower  limit  of  this 
circle  has  been  reached,  follows.  It  is 
identical  to  the  one  described  for  circle  1 
in  part  (c)  above.  All  values  of  these 
coordinates  are  assigned  to  the  appropriate 
Y(l)  and  X(I , 2) . 

(The  above  operations  are  repeated  in  the 
program  many  times.  Table  C-l  )  which  shows 
a  typical  vector  Y(I)  and  its  associated  matrix 
X(I,J)  will  be  discussed  below). 

e.  Figure  C-3  shows  two  circle  points  on  the 

horizontal  line  through  Y(I).  These  are  X(I,1) 
on  circle  1  and  X(I,2)  on  circle  2.  Since 
X(I,2)  is  closest  to  the  tooth  center  line,  it 
becomes  the  accepted  value  for  the  x-coordinate 
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X(I,J) 

st  value  in  any 


row.  All  values  in  inches). 


C-13 


of  the  tooth  profile  associated  with  Y(I). 

(In  the  actual  program  there  are  many  X(I,J) 
values  for  a  given  Y(I).  Again,  that 
x-coordinate  which  is  closest  to  the  tooth 
center  line  becomes  the  desired  point  on  the 
root  profile) . 

Table  C-l  illustrates  the  relationship  between  the 
vector  Y(I)  and  the  matrix  X(I,J)  with  selected  numbers  from  an 
actual  computer  run.  Column  1  contains  the  counter  I  which 
identifies  the  y-coordinates  on  the  fillet  circle  segments. 

While  I  takes  values  from  1  to  459,  in  steps  of  one,  in  the 
actual  run,  the  table  shows  only  a  limited  number  of  y-coordi¬ 
nates.  Column  2  contains  these  coordinates.  Each  one  is  common 
to  all  circle  points  on  the  same  horizontal  line.  Columns  3  to  8 
list  the  x-coordinates  of  these  circle  points.  With  ^  =  0°, 

i|/^n  =  -15°  and  Atp  =  -3°,  the  counter  J,  which  identifies  the 
circle  number,  ranges  from  1  to  6.  The  first  number  in  each  column 
represents  the  x-coordinate  of  the  effective  point  of  the  rack  in 
the  given  position.  (This  was  the  basis  of  the  choice  of  the  I's 
in  the  table). 

Once  such  a  "table"  has  been  completed  by  the  computer, 
the  x-coordinates  of  the  root  profile  are  determined  by  the  largest 
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(i.e.  least  negative)  number  in  any  one  row  of  the  matrix 
X(I,J).  For  example  for  I  =  231,  where  y  =  3.398  in.,  the 
profile  coordinate  x  =  -.968  in. 

3_. _ Analytical  Method  of  Obtain  Root  Profile  Coordinates 

Figure  C- 4  shows  a  rack  with  a  fillet  radius  rc  which 

is  generating  the  indicated  root  profile.  Since  the  Euler-Savary 
(**) 

Curvature  Theory v  shows  that,  the  contact  point  between  a  gener¬ 
ating  curve  and  its  envelope  lies  on  a  line  connecting  the  instant 
center  P  and  the  center  of  curvature  C  of  the  generating  curve, 
point  S  represents  a  point  on  the  root  profile.  To  obtain  the 
coordinates  of  this  envelope  point  the  components  of  the  vector 

0S  =  OC  +  CS  =  O£  +  rn  (CIO) 

c  cs 

must  be  determined.  The  unit  vector  n  has  the  direction  of 

cs 

line  PC,  as  indicated  above. 

The  components  of  the  vector  OC  were  given  by  equ's. 

(C3)  and  (C4)  or  by  equ's.  (C5)  and  (C6)  for  undercut  and  non¬ 
undercut  teeth,  respectively. 

In  order  to  obtain  a  single  set  of  expressions,  certain 
angles,  found  in  the  above  equations,  and  which  were  first  introduced 

Rosenauer  and  A.  H.  Willis:  Kinematics  of  Mechanisms,  pp.  60-62, 
Dover  Publications,  Inc.,  No.  S1796,  New  York,  1967. 
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Fillet  Radius  CF  =  rc 

Tooth  Centerline 


Figure  C-4 

Determination  of  Foot 
Profile  Coordinates  by 
Curvature  Theory 


in  Appendix  B,  are  now  redefined  as  follows  (see  Figure 
C-A): 

For  undercut  teeth,  let 


(Cll) 


[See  also  unit  vectors  R  „  and  n _  for  undercut  teeth,  as  given 

op  pr  ° 

by  equ's.  (B5)  and  (B6)]. 


For  non-undercut  teeth  ,  let 


t  =  n  = 


T 


+  0  + 


ef  f 


kjjSinO 


tan0 


(C12) 


[Compare  the  above  to  the  angles  of  the  unit  vectors  of  equ's. 
(B13)  and  (B1A) ] . 

With  the  above,  equ's.  (C3)  to  (C6)  can  be  generalized  to  the 
following  form. 

For  the  x-component  of  vector 
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(C13) 


x£  =  -  (Rp  -  b2)  sin  (t  +  \p)  +  (Rpi|;  +  a 2)  cos  (t  +  \p) 

For  the  y- component  of  vector  E5C: 

Yc  =  (Rp  "  b2)  cos  (t  +  <Ji)  +  (Rpii>  +  a2)  sin  (t  + 

(C14) 

To  derive  an  expression  for  the  unit  vector  ncs>  it  is  first 
necessary  to  express  the  unit  vector  nop  and  n  r .  alluded  to 
earlier,  in  terms  of  the  angle  t .  Thus, 

"op  '  *sin  (t  +  (',>i  4'  cos  (T  +  *)j  (C15) 

and 

npr  =  cos  (t  +  ip)i  +  sin  (t  +  ij>)j  (C16) 

Then,  as  Figure  C-A  shows: 
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n 


cs 


TC 

im 


(PR)n 


EL 


(RC)n 

°£ 


/  (PR) 2  +  (RC) 2 


(C17) 


or,  with  the  appropriate  substitution: 


n. 


cs 


(R  t|/ +  a2)[cos<T  +  t}))I  +  sin(r  +  3  -  b2[-sin(x  +  i|>)i  +  cos(x  +  4/)j] 


__E 


/  (R  ij;  +  a9)  2  +  b9 


(C18) 


or 


n. 


cs 


+  32)cos(t  +  ip)  +  b2sin(i  +  ^)]I  +  [(Rpi|i  +  a2)sin(t  +  ^)-b2cos  ( x + ^  )33 


/  a2)2  +  b2 


(C19) 


Finally,  the  root  envelope  coordinates  xenv  and  yenv  may  be 
determined  according  to  equ.  (CIO). 
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US 


(C20) 


=  xi+yj  +  rn 
c  -'c  c  cs 


with  the  help  of  equ's.  (C13)  ,  (C14)  and  (C19) ,  one  obtains 


xenv  =  " ^Rp  '  b2^  sin  (T  +  ^  +  (Rp^  +  a2^  cos  (T  +  ^ 


(R  <P  +  82)  cos  (t  +  ip)  +  bj  sin  (t  +  ip) 
r  — ° . .  ■  -  -  —  — - 

^  (Rp'i1  +  a2) 2  4-  b’ 


(C21) 


and 


env 


=  (R  -  b„)  cos  (x  +  ip)  +  (R  ip  +  a9)  sin  (x  +  <p) 
P  ^  P  *• 


+  r 


(Rp<p  +  82)  sin  (x  +  ip)  -  b2  cos  (t  +  ip) 

'  (Rp'P  +  a2)2  +  b* 


(C22) 
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APPENDIX  D 


INITIAL  AND  FINAL  ANGLES  FOR  TROCHOID  AND  INVOLUTE  GENERATION 


When  computing  gear  tooth  profile  coordinates  it  is 
necessary  to  know  the  initial  and  final  values  of  the  trochoid 
generating  angle  (see  Appendix  B)  as  well  as  of  the  involute 
generating  angle  a  (called  roll  angle  in  Appendix  A). 

The  following  shows  how  the  relationship  of  the  trochoidal 
path  of  the  effective  point  of  the  rack  cutter  with  respect  to 
the  rest  of  the  gear  tooth  profile  may  be  utilized  to  determine 
these  angles  for  undercut  and  non-undercut  teeth.  It  must  be  kept 
in  mind  that,  while  the  effective  ,'M.nt  of  the  cutter  tooth  always 
determii  s  the  magnitude  of  the  inner  form  radius  and  with  that 
the  initial  trochoid  generating  angle  ^in.  it  does  not  furnish 
the  angle  associated  with  the  lowest  point  of  the  gear  root 
whenever  there  is  a  fillet  radius  on  the  cutter.  This  final 
trochoid  angle  is  defined  as  that  angle  for  which  the  effective 

point  of  the  cutter  coincides  with  the  centerline  of  its  tro- 
choidal  path  (see  sections  4  and  5  below) . 

The  initial  involute  generating  angle  a^n  is  a  function 
of  the  inner  form  radius.  The  final  involute  generating  angle 
Ufin  depends  on  the  outside  radius  of  the  gear  tooth. 

Sections  1  and  2,  on  the  trochoid  angle  t  and  on  the 
determination  of  the  inner  form  radius  for  undercut  teeth 
respectively,  are  taken  with  minor  nomenclature  changes  from  the 
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work  of  R.A.  Shaffer:  An  Analysis  of  the  Undercut  Problem  in 
Involute  Spur  Gearing,  Report  No. R-1606,OMS5530. 11. 557 ,  DA 
Project  505-01-003,  Frankford  Arsenal,  Pitman-Dunn  Laboratories, 
Hay  1962. 


The  Trochoid  Angle  x 


Figure  D-l  shows  the  trochoid  described  by  the  effective 
point  of  the  indicated  rack  tooth.  The  line  S'U',  of  length  bg^£  , 
coincides  with  the  center  line  SU  of  the  trochoid  when  the 
generating  angle  ip '  ,  as  used  for  this  purpose,  is  zero  together 
with  the  length  QS  ,  along  the  rack  pitch  line.  The  trochoid 
angle  x  is  measured  from  the  trochoid  center  line  and  it  defines 
the  position  of  the  radius  R  (used  later  to  get  the  form  radius 
RjO  from  the  orgin  0  to  the  instantaneous  location  of  the  effective 
rack  point  at  U' .  Then, 

T  =  0  -  ip '  (Dl) 


where 


(D2) 


(D3) 
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MICROCOPY  Rl  SOLUTION  II  SI  CHARI 


Substitution  of  equ's  (D2)  and  (D3)  into  equ.  (Dl)  furnishes: 


2.  Determination  of  Inner  Form  Radius  for  Undercut  Teeth 

Figure  D-2  shows  the  side  of  a  basic  rack  of  pressure 

angle  0  in  line  with  the  base  radius  OV  =  R^,  of  the  indicated 

involute  surface.  Point  V  represents  the  theoretical  origin  of 

this  involute.  The  effective  point  of  the  rack  cutter,  shown  at 

point  U' ,  is  positioned  such  that  the  involute  profile  will  be 

undercut,  i.e.  b  >  R  sin20.  The  trochoid  of  this  effective 
er  t  p 

point  intersects  the  involute  at  point  P  and  the  length  OP 
represents  the  inne^  form  radius  R£.  Note  that  the  center  line 
of  the  trochoid  is  represented  by  line  OU,  which  forms  angle  6 


with  the  side  of  the  rack  tooth. 

To  obtain  an  expression  for  R^  one  sets  the  following 
angles  equal  to  each  other: 


t  =  6  +  8 


where 


t  =  obtained  from  equ.  (D4) ,  with  R  =  R^ 


Further , 


6  =  6  -  i|>'  . 
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(D5) 


(D6) 


where  ip’  may  now  be  expressed  in  terms  of  the  rack  pressure 
angle  6.  [  See  also  equ.  (D3) ]. 


r 


(Rp  -  bef f ) tan0 

% 


(D7) 


The  angle  6  has  its  origin  in  involutometry .  If  the  angle  a 
represents  the  involute  generating  angle  with  respect  to  base 
radius  OV,  and  if  one  defines  the  angle  POS  =  <J>,  then: 


3  =  INV  4>  =  tamp 


tan 


(D8) 


Substitution  of  equ's.  (D4) ,  (D6) ,  (D7)  and  (D8)  into  ecu.  (D5) 
furnishes : 


tan 


-i^F 


beff) 2 


eff J 


R_  -  b 


(R„  - 


=6  - 


beff)tan0 


R_ 


-  tan 


(D9) 


This  expression  may  be  solved  for  R£  either  by  a 
computer  adaptation  of  the  Newton  -  Raphson  method  or  by  a 
trial  and  error  method. 
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3. 


Determination  of  Inner  From  Radius  for  Non-undercut  Teeth 


Figure  D-3  shows  the  effective  point  of  a  rack  cutter 
making  initial  contact  with  the  prospective  involute  surface 
above  the  origin  point  V  of  the  involute.  Since  now  b^^  < 
RpSin20,  the  tooth  will  not  be  undercut. 

As  in  Figure  B-4  of  Appendix  B,  the  contact  point  Q 
and  with  that  the  effective  point  lie  on  the  line  of  action 
UP.  The  inner  form  radius  R^  may  be  described  as  follows: 


(qu)2  +  r2 


(DIO) 


where 


qu  =  PU  - 


=  R,tan0  -  ef f 
0  sin0 


(Dll) 


with  the  above,  equ.  (DIO)  becomes: 


Rc  =  /(Rutan0  -  eff  )2  +  R2 

r  D  sin0  D 


(D12) 


For  the  special  case  when  bgf£  =  R  sin20,  equ.  (D12)  becomes: 
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Rack  Side 
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Trochoid  Center  Line 


(D13) 


=  R. 


as  is  to  be  expected. 


4_. _ Determination  of  Initial  and  Final  Trochoid  Generating  Angles 

for  Undercut  Teeth 


Figure  D-4  shows  the  outline  of  the  rack  cutter  in  its 
standard  position,  as  defined  in  Appendix  B.  The  effective  point 
of  the  rack  makes  contact  with  the  undercut  trochoid,  generated  by 
it,  at  point  S.  The  flank  of  the  rack  is  in  lire  with  base  radius 
OV  =  R^,  where  V  represents  the  origin  of  the  involute  curve.  The 
associated  line  OF  =  R^  to  the  pitch  point  is  at  the  pressure  angle 
0  of  the  cutter  with  respect  to  line  T5V.  The  angle  is  zero. 

Also,  the  distance  FR,  along  the  pitch  line  and  normal  to  be£p 
eouals  a,  as  given  by  equ.  (B2)  in  Appendix  B. 

In  order  for  the  effective  point  to  make  contact  with 
point  S^n  at  the  trochoid  center  line,  the  rack  pitch  line  must 
roll  clockwise  through  the  angle  i)^n  until  points  Pf^n  and  Rfin 
coincide.  The  distance  rolled  through  equals  a  =  FF.and  therefore: 

^f in  =  '  (D14) 
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(Note  that  the  above  angle  is  negative  because  of  the  zero 
position  of  i(j.) 

When  the  effective  point  of  the  rack  cutter  is  located 

at  point  S.  ,  the  intersection  of  the  trochoid  and  the  involute 
in 

surface,  ij)  =  ij^n  .  To  determine  this  angle,  one  must  first  find 
the  total  roll  angle  Ytot  of  the  pitch  line,  associated  with  the 
change  of  the  effective  point  from  S^n  to  S^n  ,  v  depends  on 
the  arc  length  =  Pin^in  anc*  ^  equa^s  t^ie  between 

the  pitch  radii  0P^n  and  OPf^n-  With  the  inner  form  radius 
known,  one  obtains 


P.  P,. 
in  fin 


=  TS. 


(D15) 


Rj  must  be  computed  with  the  help  of  equ.  (D9) . 


Finally,  as  figure  D-4  indicates: 


or 


<J). 

*in 


/R£  -  <Rp  -  be£f)2 


(D18) 


5.  Initial  and  Final  Trochoid  Generating  Angles  for  Non-Undercut 

Teeth 


Figure  D-5  shows  the  effective  point  of  the  rack  cutter 
making  initial  contact  with  the  involute  profile  at  point  S.  Since 
this  point  lies  on  the  line  of  action  TJF,  located  above  the  origin 
V  of  the  involute,  there  will  be  no  undercutting.  (Note  that  this 
figure  describes  the  same  relationship  between  rack  and  involute 
as  Figure  B-4  of  Appendix  B  ).  The  distance  PR  along  the  rack  pitch 
line  equals  a,  as  defined  by  equ.  'BIO).  The  generating  angle 
equals  zero  in  this  initial  position,  while  the  angle  between  the 
pitch  radius  UP  =  (with  point  0  not  shown)  and  the  base  radius 
OU  =  equals  0,  the  pressure  angle  of  the  rack. 

In  order  for  the  effective  point  of  the  cutter  to  make 
contact  with  point  S^n  at  the  trochoid  center  line,  the  rack  pitch 
line  must,  as  in  the  last  section,  roll  clockwise  through  the  angle 
^fin  unt*-l  P°i-nts  P^n  and  concide.  The  distance  rolled 

through  is  again  given  by  PR  *  a.  (It  is  important  to  remember 
that  the  length  a  is  now  given  by  equ.  (BIO) ,  rather  than  by  equ. 
(B2)  of  Appendix  B) .  Therefore, 
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Ip 


fin 


a 


(D19) 


It  can  be  shown  that,  between  points  S  and  S-.  the 
trochoid  profile  moves  away  from  the  theoretical  involute  profile  SV.  In 
addition,  point  S  also  respresents  the  last  point  of  tangency 
between  the  trochoid  and  the  involute,  as  angle  ip  increases.  To 
make  this  clear,  figure  D-5  shows  that  point  S' ,  which  is 
associated  with  the  arbitrary  positive  angle  ip'  ,  is  at  a  consider¬ 
able  distance  from  the  involute  profile.  For  the  above  reason 


ip . 
m 


0 


(D20) 


6.  Initial  and  Final  Involute  Generating  Angles  for  Undercut  and 

Non-Undercut  Teeth 


Figure  D-6  indicates  that,  once  the  inner  form  radius 
is  known,  the  initial  involute  generating  angle  cun  may  be  computed, 
regardless  of  whether  the  tooth  is  undercut  or  not,  with  the  help 
of  the  following: 


ain 


T.  U. 
in  in 


(D21) 
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V 


Point  V  is 
origin  of  involute 
at  base  circle 


Figure  D-6 

DETERMINATION  OF  INITIAL 
AND  FINAL  INVOLUTE 
GENERATING  ANGLES 


Since 


T.  U. 
in  m 


*b 


(D22) 


the  above  becomes: 


a 


in 


/RI  -  g 

Rb 


(D23) 


Similarly,  with  the  outer  form  radius  R^  known,  the 
final  involute  generating  angle  is  given  by  : 


a 


fin 


Ar 


fo 


Rv 


(D24) 
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APPENDIX  E 


GEAR  STRENGTH  CALCULATIONS  .  DETERMINATION  OF  LEWIS  AND  AGMA  FACTORS 

This  appendix  first  gives  a  review  of  the  origins  of  the 
Lewis  gear  strength  formula  and  the  associated  Lewis  and  AGMA  Geome¬ 
try  Factors. +  Subsequently,  the  essentials  of  a  computer  determina¬ 
tion  of  the  Lewis  Factor  for  any  involute  gear  profile  are  developed 


For  more  detail  see:  AGMA  Information  Sheet.  Geometry  Factors  for 
Determining  the  Strength  of  the  Spur,  Helical .Herringborne  and  Bevel 
Gear  Teeth.  AGMA  226.01,  American  Gear  Manuf.  Assoc.,  Aug.  1970. 
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1.  Maximum  Bending  Stress  in  a  Parabolically  Shaped  Beam 

Figure  E-l  shows  a  cantilever  beam  with  a  parabolic 

cross  section.  The  parabola  originates  at  point  0  and  is  described 
by 

y  =  B  x2,  (El) 


where  B  represents  a  constant.  The  beam  of  thickness  b  experiences 
the  indicated  endload  F  .  Using  the  usual  bending  stress  formula, 
the  maximum  bending  stress  omax  at  the  abitrary  section  f-f  becomes: 


a 

max 


M  c 
I 


1.5 


F  y 
s  } 

b  x^ 


(E2) 


where  for  the  present  case: 

M  =  Fgy  ,  the  moment  at  section  f-f 

c  =  x  ,  the  distance  to  the  extreme  fiber  of  the  beam 


I 


b(2x)3 

12 


the  moment  of  inertia  of  the  arbitrary  cross  section 
of  thickness  b. 


If  one  now  substitutes  the  y  of  equ.  (El)  into  equ.  (E2) , 
the  following  expression  for  the  maximum  bending  stress  at  any  point 
along  the  beam  results: 


o 

max 


F 

s 


1.5  B 
b 


(E3) 
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The  above  shows  that,  the  maximum  bending  stress  in  a 


parabolically  shaped  beam  is  constant  along  the  beam. 

2_. _ Derivation  of  Lewis  Gear  Strength  Formula 

Figure  E-2  represents  a  typical  gear  tooth  which  is  acted 
on  by  a  contact  load  F.  The  line  of  action  of  this  force  passes 
through  point  L,  the  earliest  possible  contact  with  a  mating 
tooth,  at  the  end  of  the  involute  profile.  If  there  is  no  load 
sharing  with  a  second  set  of  mating  teeth,  i.e.  the  contact  ratio 
equals  unity,  this  condition  also  produces  the  largest  possible 
bending  stresses  at  any  cross  section  below  point  V,  where  the  line 
of  action  intersects  the  tooth  centerline. 

The  load  angle  6  is  defined  by  the  line  of  action  of  force 
F  and  the  normal  to  the  tooth  centerline.  The  bending  load  Fg  is 
obtained  from 

Fg  =  F  cos  8,  (E4) 

and  it  acts  normal  to  the  centerline  at  point  V. 

Now,  let  point  V  become  the  origin  of  a  parabola  which 
is  tangent  to  the  tooth  profile  at  some  section  e-e,  as  indicated 
in  Figure  E-2.  Since,  as  has  been  shown  in  the  previous  section, 
the  bending  stress  along  a  parabolically  shaped  beam  is  constant, 
section  e-e  must  represent  the  location  of  the  maximum  bending 


Figure  E-2 

Background  to  Lewis  Formula 

(Line  of  action  of  load  F,  as 

shown,  corresponds  to  unity  contact  ratio). 


stress  of  the  gear  tooth.  All  other  cross  sections  have  portions 
outside  the  parabola  and  will  therefore  experience  lower  bending 


stresses  due  to  the  load  F  . 

S 


Thus,  if  one  can  determine  the  point  of  tangency 
between  the  parabola  and  the  tooth  profile  the  maximum  bending 
stress  may  be  found  with  the  help  of  the  nomenclature  introduced 
in  Figure  E-2: 


(Fsh)(  f  ) 


6  F  h 
s 


where 

h  =  VA  ,  the  distance  from  point  V  to  the  tangency  section  e-e 
t  =  2 (AC),  the  tooth  thickness  at  section  e-e 

The  moment  arm  h  may  be  expresed  in  terms  of  the  distance 
x  =  AS  ,  which  generally  is  found  by  graphical  means.  (See  the  afore¬ 
mentioned  AGMA  Standard  226.01).  Thus,  with  the  right  triangle 
VCB  and  the  resulting  similar  triangles  AVC  and  ACB,  one  obtains  the 
relationship : 


E-6 

-116- 


and  consequently, 


h 


4  x 


(E7) 


Equ.  (E7) 
following 
tooth: 


is  now  substituted  into  equ.  (E5) .  This  results  in  the 
expression  for  the  maximum  bending  stress  in  the  gear 


o 

max 


(E8) 


In  order  to  make  equ.  (E8)  as  general  as  possible,  one 
assumes  first  that  the  dimension  x  was  determined  for  a  tooth  of 
arbitrary  diametral  pitch  P^.  The  magnitude  x^  for  a  tooth  of 
unity  diametral  pitch  will  be  given  by: 


X1  =  x  Pd 


(E9) 


The  maximum  bending  stress  in  a  tooth  of  unity  diametral  pitch 
then  becomes  according  to  the  above: 
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F 

^maxi  =  — § —  (E10) 

1  b  Y 


where  the  Lewis  factor  Y  is  given  by: 


Y 


2  x  Pd 
3 


(Ell) 


This  factor  may  either  be  directly  determined  from  a  unity 
diametral  pitch  tooth,  or  by  the  indicated  multiplication  by 
if  x  originates  from  a  smaller  tooth. 

Naturally,  for  a  given  load  F  and  tooth  thichness  b, 
^  s 

the  bending  stress  will  be  higher  for  teeth  of  larger  diametral 
pitch.  Therefore,  the  general  form  of  the  Lewis  Gear  Strength 
Formula  becomes : 


o 

max 


F  P 


b  Y 


(E12) 
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3_. _ Determination  of  Origin  V  of  Lewis  Parabola 

Before  the  computer  procedure  for  the  determination 
of  the  point  of  tangency  between  the  Lewis  parabola  and  the  tooth 
profile,  as  well  as  of  the  distance  x  ,  can  be  given,  it  is  first 
necesary  to  find  the  distance  Rv  =  0V  from  the  center  of  the  gear 
to  point  V,  the  origin  of  the  Lewis  parabola. 

Figure  E-3  shows  the  line  of  action  of  force  F  passing 

through  the  arbitrary  point  WR  of  the  involute  profile.  (Since  the 
force  F  is  always  normal  to  the  involute  profile,  its  line  of  action 
coincides  with  line  URWR  ).  The  radius  R  =  0WR  is  associated  with 
point  WR  . 

The  load  angle  BR  of  this  line  of  action  with  the  normal  to 
the  tooth  centerline  is  identical  to  the  angle  between  the  tooth 
centerline  and  the  indicated  base  radius  0UR  =  R^  .  This  later  angle 
may  be  defined  as  follows: 


eR  *R  "  0R 


(E13) 


where 


(E14) 
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I 


and 


(E15) 


The  above  arc  length  is  equal  to  one  half  of  the  circular 
tooth  thickness  at  radius  R.  To  determine  this  tooth  thickness, 
one  makes  use  of  the  following  relationship  from  involutometry 
which  allows  the  determination  of  a  circular  tooth  thickness  T2  at 
a  radius  R2  if  the  tooth  thickness  T^  at  radius  is  known: 


T2  "  Tl-ld 


-  2R2(  INV  <j>2  -  INV  ^  ) 


(E16) 


where 

INV  <t>i  =  tan  <t>i  -  <i>i  (E17) 

To  adapt  the  above  to  the  present  case,  it  is  necessary 

to  know  the  actual  circular  tooth  thickness  T  at  the  standard 

cs 

pitch  radius  Rp  of  the  gear.  (The  magnitude  of  Tcg  depends  on 
whether  the  gear  is  standard  or  modified) .  The  associated  angle 
4  corresponds  always  to  the  pressure  angle  0  of  the  rack  cutter. 
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For  the  case  when  the  contact  ratio  of  the  mesh  is  unity 
and  therefore  the  crucial  contact  is  made  at  the  end  of  the  involute, 
one  lets  R  =  Rq  ,  the  outside  radius  of  the  gear,  in  eau.  (E20) .  This 
gives  the  load  angle  6  for  the  Lewis  parabola: 
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R.  T 

6  ■  C°S'I(  r  >  '  5g  +  INV  V 


inv  e 


(E21) 


where 


-1 

cos 


(E22) 


To  determine  the  position  of  point  VR  ,  where  the  line 
of  action  intersects  the  tooth  centerline , one  considers  the  following 
(See  Figure  E-3) 

\  *  WR  *  «R  -  W  <E23) 


Since  the  line  WCM_  is  normal  to  the  tooth  centerline, 


WrMr  =  Rsin0R  , 


(E2A) 


while 

0Mr  =  Rcos0R  (E25) 
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(E26) 


MrVr  =  R  sinOR  langR 


Substitution  of  equ's.  (E25)  and  (E26)  into  equ.  (E23) 
leads  to  the  position  of  point  VR  with  respect  to  center  0  of  the 
gear: 


Rv  =  R  (  cos -  sin0Rtan8R) 


(E27) 


For  the  case  when  R  =  Rq ,  the  above  expression  becomes 


R  =  R  (cosQd 
v  o  R 

o 


sin0„  tang) 
Ro 


(E28) 


where  B  is  given  by  equ.  (E21)  and  0R  is  obtained  with  the  help 

*0 

of  equ.  (E19)  ,  i . e . : 


cs 

2Rn 


INV  <*>„  +  INV  0 

Ro 


again,  <J>R 


-1  Rb 
=  cos  w- 

Ro 


(E29) 


(E30) 


Computer  Determination  of  Lewis  and  AGMA  Geometry  Factors 


The  computer  determination  of  the  Lewis  and  the  AGMA 
Geometry  factors  requires  that  the  coordinates  of  the  gear  tooth 
profile  are  fully  known  and  stored  in  the  associated  vectors 
X2(K)  and  Y2(K).  Subsequent  to  this,  the  following  operations 
are  performed  (See  Figure  E-4) : 


a.  The  distance  Rv  =  OV,  of  the  origin  V  of  the 
Lewis  parabola,  is  found  with  the  help  of  equ. 
(E28) . 


b.  To  accommodate  the  specific  situation,  the 
equation  of  this  parabola  becomes: 


\  -  Y2(K) 


(E31) 


where  X^ar  is  the  x- coordinate  of  the  parabola, 
associated  with  the  coordinate  Y2(K).  Again,  as 
in  equ.  (El) ,  B  represents  the  constant  of  a 
specific  parabola.  Since  Xpar  is  required,  ecu. 
(E31)  becomes: 


Y2(K) 


(E32) 
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c.  In  order  to  find  the  tangent  parabola  and  with 
that  the  value  of  the  parameter  x  of  the  Lewis 
formula,  the  appropriate  constant  B  must  be 
determined.  To  this  end,  a  comparatively  large 
B  is  first  chosen  to  make  sure  that  the  associated 
parabola  does  not  contact  the  tooth  profile.  The 
inner  parabola  of  Figure  E-4,  with  the  constant 
B^,  represents  such  a  case.  As  the  magnitude  of 
B  is  progressively  decreased,  the  resulting  para¬ 
bolae  approach  the  tooth  profile  more  closely. 

(The  constant  B2  of  the  outer  parabola  of  Figure 
E-4  is  smaller  than  the  constant  B^  of  the  inner 
one).  Eventually,  the  tangent  parabola  will  have 

the  correct  value  of  B. 


d.  The  above  process  is  accomplished  by  computing 


for  a  given  value  of  B  all  X  as  functions  of  the 
0  par 


Each  Xpar  is  compared  with 


full  range  of  the  Y2(K) 
the  corresponding  value  of  X2(K)  of  the  tooth  profile. 
Whenever 


X2(K) 


X  | 
par  ' 


(E33) 


the  correct  value  of  B  was  used  and  the  associated 
values  of  X2(Kfcan)  and  Y2(Ktan)  lie  on  the  tangencv 
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cross  section  of  the  gear  tooth.  The  tooth 
thickness  t  as  defined  by  Figure  E-2  then 
becomes : 


t  2  [ |  X2(Ktan)  |]  (E34) 


e.  Finally,  the  Lewis  distance  x  of  equ.  (E6)  may 
be  found  by: 


x  - 

h 


X2(Ktan)] 


Y2<Ktan> 


(E35) 


The  AGMA  Geometry  factor  obtained  from 

the  substitution  of  the  above  into 


YAGMA 


_ 1 

cos  g  /■! .  3 

COS0  x 


tanB  _ 

2[|X2(Ktan)|] 


(E36) 
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APPENDIX  F 


PROGRAM  MATRIX 

The  present  appendix  describes  and  lists  program  MATRIX 

which  is  used  for  the  determination  of  the  sizes  of  the  various 

arrays  in  programs  LEWIS  CIRCLES  and  LEWIS  ENVELOPE.  A  sample  of 

the  program,  together  with  two  output  pages  is  given  in  section  5 
below . 

1 .  Determination  of  Length  of  Vectors  XI  and  Y 1 

(Used  in  programs  LEWIS  CIRCLES  and  LEWIS  ENVELOPE . ) 

The  vectors  XI  and  Y1  are  used  for  the  initial  storage 
of  the  involute  coordinates  of  the  tooth  profile.  The  common  length 
of  these  vectors  is  determined  with  help  of  equ.  (1)  i.e. 

Yl-rows  =  af  in _ ai  n  (FI) 

Aa 

The  explanation  for  the  above  nomenclature  follow's  equ.  (1)  in 
section  2. 

2 .  Determination  of  J3 iz es_ o  f_ Vector  Y  an d  Array  X 
(Used  mostly  in  program  LEWIS  CIRCLE.S.) 

The  number  of  rows,  both  of  the  vector  Y  and  the  array 
X,  depends  on  the  number  of  one  one- thousandth  of  an  inch  contained 
in  the  vertical  distance  between  the  initial  point  of  the  first  fillet 


F- 1 
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circle  and  the  final  point  of  the  last  one. 


The  y-coordinate  y.  of  the  first  point  of  the  first 
J  in 

circle  is  that  of  its  effective  point.  It  is  determined  for  \p  =  tp^ 
according  to  equ.  (B9)  or  equ.  (B25): 

yin  =  (Rp  -  beff)cos(T+  ^in)  +  (Rp^in  +  a)sin(t+  ^in)  (F2) 


where,  for  an  undercut  tooth 

t  =  n  =  J'1  1  0,  according  to  equ.  (B9)  (F3) 

I 

y^n  =  initial  roll  angle  of  rack  pitch  line  according  to 
equ.  (D18) 

For  a  non-undercut  tooth: 


T 


n 


+  0  + 


b 

R 


eff 

p  sinG 


tanO,  according  to  equ.(B?5)  (F4) 


=  0,  according  to  equ.  (D20) 

The  final  point  y£^n  of  the  last  circle  is  obtained 
according  to  equ.  (3)  or  (5),  i.e.: 
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-^fin  ^cfin  ”  rc  (F5) 

ycfin  represents  the  center  of  the  fillet  circle  for  \ (i  =  ^fin- 
Then,  according  to  equ's.(C4)  or  (C6) : 

yCfin  =  (Rp  "  b2)cos(t+  ^f.n)  +  CRp^fin+  a2)sin(x+  *fin) 

(F6) 

When  the  tooth  is  undercut  the  angle  x  is  used 
according  to  equ.  (F3) ,  otherwise  equ.  (F4)  must  be  used.  The 
angle  i|^n  for  undercut  and  non-undercut  teeth  is  the  same.  It 
is  given  by  equ.  (D14)  or  equ.  (D19) . 

Thus,  the  number  of  rows  of  X  and  Y  are  found  from: 

Y-rov/s  =  1000(y^n  -  Yf^n)  [in  form  of  an  integer]  (F7) 

The  number  of  columns  of  the  array  X  depends  on  the 

number  of  circles  involved.  This  in  turn  is  a  function  of  the 

angular  increment  A vp  as  well  as  and  . 

°  in  fin 

Therefore,  as  in  equ.  (10): 

X-columns  =  ^  ^  [also  made  into  an 

integer] 

3.  Hand  Computation  of  Sizes  of  Vectors  X2  and  Y2 

(Used  in  programs  Ll'Wfs  CIRCLF.S  and  Ll-v  IS  F.NVFI.OPV) 
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The  size  of  vectors  X2  and  Y2,  which  contain 
the  final  tooth  profile  coordinates,  are  not  computed  in 
program  MATRIX,  but  may  be  obtained  by  adding  the  results  of 
equ's.  (FI)  and  (F7) : 

Total  rows  of  X2  and  Y2  =  Yl-rows  +  Y-rows  (F9) 


The  above  will  give  an  exact  result  for  use  in 
program  LF.KTS  ClPCLF.S  and  a  somewhat  larger  than  needed  number 
for  program  LFWIS  FNYll.OPF.  (The  y- coordinates  of  the  root  are  not 
generally  computed  every  one  one- thousandth  of  an  inch  hv  this  routine.) 


4.  Description  of  Computer  Program  MATPIX  (See  section  5  of  this 
Appendix) 


a.  Input  Parameters 

The  following  input  data  are  required: 


N 


RIFF 


PD 

TUFTAD 

CAPTCS 


number  of  teeth  in  gear  under  consideration 
beff  for  unity  diametral  pitch  (See  Appendix  B) 

=  1.000,  the  diametral  pitch 
9,  the  pressure  angle  of  the  hob 
Tcg ,  the  circular  tooth  thickness  of  the  gear  at 
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RC 

KADD 


DHL PS  I 


DELAL 


the  standard  pitch  radius 


•in  terms  of  Pj  =  1.000 


r  ,  the  rack  fillet  radius  for  unity  diametral  pitch 

Addendum  constant  of  gear,  with  the  magnitude  of  the 
addendum  defined  by  KADD/ PD,  in  terms  of  Pd  =  1.000 


A\p ,  the  increment  of  the  roll  angle  of  the  rack 
pitchline  on  the  pitch  circle  of  the  gear 


Aa,  the  increment  of  the  involute  generating  angle 


b .  Computations 

The  initial  computational  sequence  in  program  MATRIX 
is  in  many  ways  identical  to  that  of  the  early  parts  of  programs 
LEWIS  (See  section  2).  Following  the  test  which  establishes 
whether  the  tooth  under  consideration  is  undercut  or  not,  the 
appropriate  parameters  for  equ's.(Fl)  to  (F8)  are  computed. 

The  evaluation  of  these  expressions  starts  in  line  61 


i  n  the 

following  sequence  : 

YIN 

=  Yin>  according  to  equ. 

(F2) 

YCFIN 

=  ycfin’  accordin8  to  e9u 

•  (F6) 

YFIN 

=  Yfin'  according  to  equ. 

(F5) 
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The  integral  number  of  rows  of  X  and  Y  is  found 
in  lines  64  to  66,  according  to  equ.  (F7).  The  integral  number 
of  columns  of  array  X  is  determined  according  to  equ.  (F8)  in 
lines  67  an  68.  The  number  of  rows  of  vectors  XI  and  Yl,  based 
on  equ.  (FI),  are  found  in  a  similar  manner  in  lines  69  and  70. 


c .  Output  Parameters 

(All  lengths  are  given  in  inches  and  all  angles  are 
given  in  degrees) 


The  output  of  the  program  first  states  the  numerical 
values  of  the  input  parameters  N,  BFFF,  PD,  THF.TAD,  CAFTCS,  RC , 
KADD,  DELPSI,  and  DFLAL . 

Subsequently,  it  is  noted  whether  the  tooth  is  undercut  cr  rot. 
In  addition,  the  numerical  values  of  the  following  computed 
quantities  are  given: 


RF 

ALFHIN 

ALFHFIN  = 

FSIIN 

PSIFIN 

YIN 

YCFIN 


in 


“fin 


ip  ■ 
rin 

ipr. 

Yf  in 

^in 

^efin 
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YFIN 


in 

Finally,  the  number  of  rows  of  the  vectors  XI  and  Y1 , 
together  with  the  size  of  the  tooth  root  array  X,  are  printed  out. 

5 .  Sample  Program  and  Sample  Output 

The  following  pages  contain  a  listing  of  program  MATRIX 
as  well  as  two  sample  outputs. 

The  programs  were  run  for  the  same  teeth  as  found  in 
programs  1  i-KTS  CIRCl.FS  and  LFWIS  ENVELOPE  in  section  3  of  this 
report . 

The  input  data  for  the  first  run  on  p.  F-12  are  as  follows 


N 

- 

36, 

BEFF 

= 

1.215  in. , 

PD 

= 

1.000, 

THF.TAD 

= 

2  0°, 

CAPTCS 

= 

1.8460  in. , 

RC 

= 

0.040  in., 

KADD 

= 

1.425  in., 

DEEPS  I 

= 

0.125  degrees 
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DliLAL 


0 . 25  degrees 


■  The  program  found  that  the  tooth  is  not  undercut.  The 
various  output  parameters  from  RV  toYFIN  may  be  used  for  checking. 
Finally  the  array  sizes  are  given  as: 

Length  of  Vectors  XI  and  Y1  94 

Length  of  Vector  Y  :  370 

Size  of  Array  X  370  x  85 

The  lengths  of  vectors  X2  and  Y2  are  obtained  by  addition 
of  the  lengths  of  Y1  and  Y,  according  to  equ.  (F9) ,  i.e. , 


94  +  370  =  464 


The  second  run  of  the  program  on  page  F -  3  3  concerns  itself 
with  an  undercut  12  tooth  pinion. 
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$4*$$  V  4  T  R  I  X  4  4  4  4  4 


1 


A 

5 

6 
1 
8 
9 

1  0 
I  1 
1  2 
1  3 
:  a 

>  c 

I  6 
l  7 
18 

1  <3 

:c 

:  i 
>2 
?  3 
?A 

25 
76? 

2  3 

26 


2S 
3  C 

31 

32 

33 
3A 
3* 
36 
3  7 
36 
3C 
AC 
A  1 


A2 
A3 
<  A 
A  5 
A  C 


t  JOB 
C 

c 

c 


IMPLICIT  KFAt*8( A-H.C-7)  i 

C  C  *>  N  C  N  RF ,P  ,RP .THETA  j 

REAL  N.KACC 

3 E AC  (5,5  1  )  N.BFFF.PD.THFTAD.CAPTGS.RC 
RE  AC ( 5 , 52  )  KACC 
RF  AC  ( 5,5?)  CFLOSI.DELAL 
51  F  OR  v  *  T ttF  3.5) 

52Fn3VJT(F3.5) 

53  FQi3*'.*T(2F9.5) 

V» R  1  T H  (6  ,  l  )N»PEFF,PD 

1  FORMA  I  (  •  1  •  ,  *\  =  ’,F9.5,r>X,’BEFP  =  •  ,F9.5,5X,  'PC  =',F9.5) 
'.>RITF(6,3)THcTAC,CAPTCS,RC 

3  FC«NM//'  THETAC  =  •  ,  F  9 . 5  »  5  X  ,  ' C  AP  TC S  =  •  ,  F  9 . 5 , 5  X  ,  •  R  C  *’,F9.5) 

VvR  I  TF  ( t  «  2  )  KACC.DELPSI  .CtLAL 

2  FORMAT (  /  /  *  '.'KACC  = • , F 9 . 5 , 5 x , • CE L P $  1  * • , F9. 5 , 5* , • DEL AL  =',F9.5//) 
oi  =  3 .  I  A l 55 

?  =  FI/1FC. 

THETA  -  T  H  E  T  AC  *  Z 
3  P  -  N/(?.*PC) 

RB  =  3F*CCCS( THETA) 

50  -  RF+':ACC/°C 
(3  -  n  E  F  F  /  7  0 

THETA6  -  CAPTCS*nCnSITHETA)/RRH2.*(DTAN(THETA l-THETA) 

FPS  =  THET  AC/2.C0  ♦  THETA 

CELIA  s  F/ (RP40S  I  N  <  THET  A  )  )  -  DTANMTHEM) 

PALLCVs  =  R  P4CS  IN { THET  A  )  *0S I M  THET  A  ) 

IFfEOFF.CT.PALLCw)  GC  TO  Al 
IFIREFF. IE. FALLOW)  GC  TC  A2 

C  44444444444444444444444444* >4"  44  4  4 

c  4444*  TCCTH  IS  UNDERCUT  444*4 

C  44444444  4  444*444444  4  444444444444 

A l  A  =  (RP-B ) 4CTAN ( THETA  ) 

TAL  =  EPS 
PT  =  RC 
C  =  RE 

RF=CEKKFR ( PT ,C , 1.0-15, I.C-IP, 100) 

ALPHINI  =  C  SORT  (  0  APS  (  RF4RF-RP4RR  )  ) /RP 
ALPHF  =  CT  AMC  ARCCSIRP/RC)  ) 

G  A  V  T  ^  (CSCRT(CARS(RF4RF-(RP-BEFF)  442  )  ))/RP 
P  S I  I N  =  -A/RP  4  GAM! 

PSIF1N  =  -A/RP 
WR  I  T  E  (  6 , 6  5  ) 

65  FORMAT!//’  44444  TCOTH  IS  UNDERCUT  44444’//) 

GO  TO  A3 

C 

C  444444444*4444 4 4 444*44 4* 44444444444*4 

c  44444  TOOTH  is  NOT  UNDERCUT  44*4*4 

C  *44**44*4*4*44*44*44****4*4444******* 

C 

A2  A  =  R/CTAN! THETA  ) 

TAU  =  EPS  ♦  CELTA 

RF  =  OSQRT(f)A0S((R0  +  DTAN(THETA).-B/DSIN(THFTA))**2*RB*iRB)) 

ALPHIN  =  OSORT (0AQS(RF*RF-RR*RB)1/R8 
IF  ( ALFHIN.LE.O.CUI )  ALPHIN  =  CT 
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ooooooooooooo  o  "i  n  n  n 


HPHF  -  CTAN(CARCOS{RR/RC>) 

F  S I  IN  =  C.O 
CSJFT\  =  — 

'<*  R  I  T  E  (  6  » t  fe  ) 

66  FCRMiTI//*  <*♦**  TCOTh  IS  NOT  UNPEROUT  ««***•//) 

43  CONTINUE 

V8ITE<6,62)  R  F 

62  FORMAT!//*  »,*RF  =*,F9.S//» 

$****•}****.«*«**>**  -Xt  <■»***«**««***«<<«>*#$ 

rr  <  /'LL  TVPES  °f  TEETH  ««*** 

<**<*<'**$V$**«***$V'$v 

f  s  1 1  n  l:  =  psiin/z 

PS  r  F  C  =  PS  IF  IN tl 
•A  L  l  N  C  =  a  L  DF  IN/? 

ALFIM'  =  ALPt-F/7 

"2  =  R  -  :'.C*CS1MTHCTA) 

A  2  =  A  —  9  C  *  C  C  0  S  ( T  H  F  T  A  ) 

YIN  -  I RO-R) *CCCS ( TAU  +  PS I  IN )+ (RP*PS I  I N ♦ 4  I *U$ INC TAu+PSI IN ) 

YCF  I  A  =  (  RP-E2  )  *nccs  (  T  AIM  PS  IFIN)4-(RP*PSIFlN+A2)«'f'S!N(  TAU+PSI FIN) 
YFIN  =  YCFIN  -  8C 
ROW  =i  Y  LA  -  YFIN 
R  C  V>  =;  ROW*  1CCG 
IRCW  =  RC’n 

COL  =  (PSIIND  -  PSrFD/OELPSI 
I CrL  =  COL 

R  D  In  I  N  V  =  (  ALPFF-ALPF  I N ) / » Z*CEL A L ) 

I  R  I  A  V  =  RCVvINV 
V>RITF(e,j)  ALINC,  ALF  INO 

5  FORMAT! //•  ‘.MLPFIN  =•  ,F9.5f5X,*ALPHFIN  =*tF9.5) 

4R  I  Tl  (6  ,86  >  PSIINC.PSIFC 

S6  FORMAT  (  /  /  *  *,*PSIIN  a*  ,F9.5,5X,  'PSIFIN  =*,F9.5) 

>  R  I  T  F ( f  , ?  2 )  YIN, YCFIN,  YF  IN 
t»RITF(E,?5)  IRINV 

P 5  FORMAT ( /////•  *  *  5X ,’ LENGTH  CF  ROWS  FOR  INVCLLTE  VECTORS  XI  AND  Y1 
l,//'  RCWS  =  *  ,  I  6  ) 
lnRITE(t»87) 

87  FORMAT (//•  ',5X,'NUvRER  CF  ROWS  AND  COLUMNS  FOR  TOOTH  RCCT  ARRAYS 
l  I 

V-RITE(6,83)  IRCW.ICCL 

fl 2  FORMAT ( //• OY l\*« ,Fll.5,5X, • YCFIN** ,Fll .5, 5X, *YFIN*»  »Fl 1.5) 

P  3  F  CRM  ATI//'  •,*  RCN*»  ,  I6.5X,  •  C0L=',L6) 

STOP 

ENC 

CCl'BLE  PRECISION  FUNCTION  DE KKE R I B , C , E PS , TCI , L I M I T ) 


TITLE-  CEKKER'S  ALGORITHM 

PURPOSE-  TC  FINC  A  RCCT  CF  A  NON-LINEAR  ALGEBRAIC  EQUATION  F (X) =0. 
DESCRIPTION  CF  PARAMETERS 

P  =  THE  LATEST  ITERATE  AND  CLOSEST  APPROX  I MA T I  ON  TO  THE  ROOT. 

A  »  THE  PREVIOUS  ITERATE. 

C  »  THE  PREVIOUS  CR  AN  CLCER  ITERATE. 

COMMENTS 

AT  ALL  TIMES  6  ANC  C  BRACKET  THE  ROOT. 

INITIALLY  A  IS  SET  EQUAL  TO  C. 

F-10 


C  LSE  EPS  =  !.D—  7  AND  TCL=1.C-10. 

C  SUPRCUTINSS  ANO/CR  FUNCTION  SUBPROGRAMS  REGUIREC-  FfX) 


ft 

f  1 
t  c 

MS 

so 
o  i 
52 
02 
9  A 

C,  e 

se 

97 
EE 
SS 
If  C 

ic  r 

If  2 
1C2 
1  C  A 
ICE 
U  fc 
!  r  7 
HE 

I  C  9 
11C 

I I  1 

112 

112 

11  A 
115 
11C 
1 1 1 
tie 

119 
1  2  C 

12  1 
122 
122 
1  2  / 

125 

126 

l  ?  1 

ire 


IMPLICIT  REAL*fi  (A-H.o-Zl 
»=C 

* CLN T  =-  l 
FCA  =  F  (  A  ) 

FO0  =  F ( E  ) 

FCC=F(C)  _ 

I F  c <FCP*FCC) .GT.O.DOIGO  TO  6 
fl  KOLNT  =  KOLNT-*'l 

IF (DAPS (FCP).LE.OABS(FCC) )GC  TO  7 
T  =  P 

e=c 

CEKKER=B 
C  =T 
A  =  C 
T  =  F  C  P 
FCn=FCC 
F  OC  =  7 
FO«=FCC 

7  IF  IFCP.EC.Q.HC 1RETURN 

T0L1=2.CC*EPS*CABS(B)+0.5D0*TCL 
IF (DATS (F-C) .LE.2»nc*TCL 11RETURN 
IF(DAPS(FCP-FC4).LT.1.C-70)RETURN 
IF(KCUNT.GE.LHIT)CC  TO  12 
n=A-{  (b-a)«f:a)/(fcr-foai 
pp^ci 

IF ( .NOT. ( ( (  C  .  L  E  .  D I 1 . INC. (Cl •  LE«  B) ).0R. I ( B. LE . 0 1 )  .  ANC.IDI.LE 
1  PPMfHCI/2.CC 

IF  C  n  (C.LF.CI  ).  AND.  I  V.  I  .  LE  .  R  ))  .OR  .  (  I  R.  L  E  .  D  I  l.ANP.  (CI.LE.C  )  )  ) 
1  . ANC. IOAPStCI-P).LT.TCLl))PP=B*nSlGN(TOLlt(C-R} ) 

IF  (CAPS  (  III-Cl.LT. CAPS!  R- Cl/A.  CO  )PP  =  (R  +  C)/2.rC 

A  =  P 

f CA=FCP 
P=PP 

CEKKERsfl 
F  OP  =  F ( F  ) 

IF((FCe*FCC).OT.O.OO)C=A 
IF( (FGHPFCC ) ,GT.0.n01FQC=F0A 
CO  TC  E 
fc  WRITE(6,15) 

£  TCP 

12  V*  R  I  T  E  C  6 , 13  » 

'TCP 

15  FURKAT ( •-» , 'FRROR  IN  INPUT  CSTA.  FIB)  HAS  TEE  S-AME  SIGN  AS 
’  '  F (C 1 • ) 

17  FORM  AT  (•-•, 'NUMBER  OF  ITERATIONS  HAS  EXCFECFC  LHIT.  STOP 
’  7  ICNS.  •  ) 

ENC 


,o  m 


C,  0  *■’  P  U  T  4 


irs 
1  2  C 
in 
1  ?  2 


1  ?  2 

1  7  A 


CCLEIE  PRECISION  FUNCTION  FIX) 

I mpl  LC 1 T  REAL*0  ( A-h, 0-Z  1 
COMCN  RF,»,RP,  THETA 

F =r A  TAN ICSCRT ( X**?- ( RP-R  )**2)/(RP-B  ) 1-DSORT I X**2- (RP-R  )**2)/RP 
1-TFETA«(MP-P  )  *CT  AN  (  THETA  l/KP-OSQRT  (T**2-RP«*2)/RP*DATA\(nSCRT(X 
2<*?-7r#->2)/RP) 


RETURN 

ENC 
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00000*1  -  Od  O9ZS0M  =  d*3d  DOiij*'! 


CORE  ISARF  CKJECT  COCE*  Bl”>7  BYTES, ARRAY  ARE  A=  4R  BYTES, TOTAL  AREA  AVAILABLE*  3B9I2  BYTES 

CIAC'.C  ST  ICS  AL  *BF  R  CF  ERRORS*  0,  NUMBER  OF  WARAlNCS*  0,  NU“BER  CF  EXTENSIONS*  0 


APPENDIX  G 


PROGRAM  LEWIS  CIRCLES 
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!AMtM  XL  ICCIt-  PKLULE 


•  <  —  -t  J-  ~ 
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<  « 
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tf  *  u-»  *-  y- 

arzanu-" 

K  1  L)  u  •  ’J  w  • 

— «  —  — .  UL  O  ^  •« 

_)  2  v/l  V#  XI  X  5fc  «»  » 

—  o  z  ^  —  * 

-**  J  u,  ,J  w  J  □  , 

a  :•.  i  i  «  i 

-uu  3  ■«  a  a  x* 


’  o  x  <  □ 

»  3  •  -J  •  C 

:  <  •  jj  *  v 
>  *  •  a  • 


:  >—  «i  x  <*  >. 

.  v-  a.  a  . 

x 

I  1  U  X  U  1 

.  -»■  u.  XX 


—  — j  x»  *■* 

XT  .1  • 

:  :j  —  >>  ‘'N  ~ 

•  *  -*»  t  a  c  i- 

i  *  T  «%  *. 

•  J  <  X 

'  \  U  t  \  _j 

-  f  *  *  fl  I  I 

N  C  I*  *.  *— 

-£  O.  -A  —  -" 

■ '  ~  'I  - 


i  U  ^  X  *v  x 


IJ  Ll  w  u  vj  o  J  w>  '  J 


_<  r\j  rj~(  .O^CO^O— 

1>‘3'30'3,JOOOOOOOOOOOO'3'^-300^*-,^no 
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fCHIKAN  IV  i  Uv£L  £1  MIA  CATE  *  81135  16/18/C4  PACE  0004 


•  |  «l  U  u  "  ^  •* 
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